#7 
MAS 


athematical Tables 


and other 


Aids to Computation 


UNIVERSITY 
OF MICHIGAN 


NOV 26 1952 


MATHEMATICS 
| 1RRARY 





A Quarterly Journal 
Edited by 


E. W. CANNON F. J. MURRAY 
C. C. CRAIG J. TODD 


A. ERDELYI D. H. LEHMER, Chairman 





VI - Number 40 - October, 1952 + p. 207-268 


Published by 


THE NATIONAL RESEARCH COUNCIL 
Washington, D.C. 











NATIONAL RESEARCH COUNCIL 
DIVISION OF MATHEMATICS 


EDITORIAL COMMITTEE 
E. W. Cannon (E.W.C.), National Bureau of Standards, Washington, D.C. 
Automatic Computing Machinery [ACM ]. 


C. C. Craic (C.C.C.), University of Michigan, Ann Arbor, Mich. Mathe- 
matical Statistics [K]. 


A. Erpétyi (A.E.), California Institute of Technology, Pasadena, Calif. 
Higher Mathematical Functions [L]. 


F. J. Murray (F.J.M.), Columbia University, New York, N. Y. Other 
Aids to Computation [OAC]. 


J. Topp (J.T.), National Bureau of Standards, Washington, D.C. Numer- 
ical Methods. 


D. H. LeEaMER (D.H.L.), Chairman, National Bureau of Standards, Los 
Angeles 24, California. 


SUBSCRIPTION RATES 
1943-1945: (Nos. 1-12) $12.00 for all 12 issues (not available for separate 
sale except as noted below*) 
1946-1949: $4.00 per year 
1950-1952: $5.00 per year 


Single issues are available for sale as follows: 


* 1944 (No. 7, ‘““A Guide to Tables of Bessel Functions,’’ by H. Bateman and 
R. C. Archibald, 104 pp.) $2.00 


1946-1949 (Nos. 13-28) $1.25 for single issue 
1950-1952 (Nos. 29-40) $1.50 for single issue 


All payments are to be made to National Academy of Sciences, 2101 Con- 
stitution Avenue, Washington, D. C. 


Agents for Great Britain and Ireland (subscription 42s, 6d for 1952) Scien- 
tific Computing Service, Ltd., 23 Bedford Square, London W.C.1. 








Published quarterly in January, April, July and October by the National Research Council, 
Prince and Lemon Sts., Lancaster, Pa., and Washington, q 


All contributions intended for publication in Mathematical Tables and Other Aids to Compu- 
tation, and all Books for review, should be addressed to D. H. Lehmer, 16556 Chattanooga 
Place, Pacific Palisades, Calif. 


Entered as second-class matter July 29, 1943, at the post office at Lancaster, Pennsylvania, 
under the Act of August 24, 1912. 





























THE ONR RELAY COMPUTER 

















The Office of Naval Research Relay 
Computer 


1. Introduction. This paper describes a general purpose relay computer. 
This digital machine was installed in Staughton Hall on the campus of The 
George Washington University on May 9, 1951. At that time it was renamed 
the ‘‘ONR Relay Computer”’ and was made available by the Office of Naval 
Research to the Logistics Research Project for operation and maintenance 
under an ONR-GWU contract. The machine is a slow-speed, single-address 
computer with magnetic drum storage of 4094 numbers or instructions, 
each consisting of 24 binary digits including sign (equivalent to about seven 
decimal digits). It is capable of performing 42 arithmetic, logical, transfer, 
output and stop instructions through use of 734 mechanical relays and 655 
electronic tubes. Input to the computer is on mechanically sensed, seven- 
level teletype paper tape, while output may be on electric typewriter, or on 
punched tape. A view of the computer is shown in the frontispiece. 

2. Storage. The storage element is a cylinder or drum 12 inches in diam- 
eter, 83 inches long. The cylindrical surface of this drum is coated with a 
ferromagnetic material, small areas of which may be magnetized in either 
of two directions to represent the binary digits 0 or 1. The drum rotates 
at 440 revolutions per minute, and the magnetization configurations can 
be played back from or recorded (stored) on it 24 at a time by a group of 
24 magnetic heads placed .002 inches from the surface. Two additional 
playback heads are required for timing. All electronic tubes in the computer 
are used in connection with the storage unit; they accomplish amplification, 
shaping and timing of signals in the memory unit. These groups of 24 binary 
digits represent both numbers and coded instructions to the computer, 
either of which may be stored on any part of the drum. Successive memory 
locations are normally loaded from punched paper tape sensed by a stand- 
ard mechanical tape reader, but any can be loaded singly by manual controls. 
Each group of 24 binary digits has an address of 12 binary digits which 
refers to and specifies its particular location on the drum. There are 4094 
such addresses, from octal 0002 through 7777. Given any address, the 
operator can inspect by manual control the contents of the address in a 
display of neon lights. 

3. Arithmetic Unit. All arithmetic is performed by mechanical relays 
within the computer. A relay may be either open (nonconducting) or closed 
(conducting), representing the binary digits 0 and 1, respectively. The 
computer contains three basic arithmetic registers which are groups of 
relays: the Accumulator, the Q register and the X register. The Accumu- 
lator, A, is capable of binary representation of a number of 48 digits. It is 
the principal arithmetic register, the place where the actual arithmetic is 
performed, and possesses subtracting and shifting properties. The 24-binary- 
digit Q register holds the multiplier during mu!tiplication and the quotient 
after division. The 24-binary-digit X register serves as a temporary storage 
for both numbers and commands as they come from the memory. In it are 
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208 THE OFFICE OF NAVAL RESEARCH RELAY COMPUTER 
held instructions during decoding and data while being operated upon. For 
example, the X register contains the multiplicand during multiplication. 

The binary number system is used throughout. On the control panel 
lights and buttons are grouped in triplets for easier octal representation. 
All numbers are considered as integers by the computer, with the left-most 
digit indicating the sign (0 designates positive; 1, negative). A negative 
number is represented by the ‘‘one’s” complement of the corresponding 
positive number. The Accumulator is basically a subtractor, so additions 
are actually performed by complementing one of the numbers to be added 
and subtracting it with end-around borrow from the other number. Such 
subtractions are performed in parallel; i.e., all 24 digits are operated on at 
once. Multiplication is a series of subtractions and shift operations. 

Conceptually, the arithmetic unit is much like a desk calculator where 
the Accumulator corresponds to the product or upper dials, with twice as 
many digits as the keyboard; the Q register corresponds to the lower or 
multiplier dials, equal in length to the keyboard; and the X register is 
similar to the keyboard itself. But the Relay Computer’s ‘keyboard’ X 
can effectively be filled from any of the 4094 addresses under control of 
instructions on the drum. Basically the arithmetic unit is like an auto- 
matically sequenced desk calculator with 4094 keyboards. 

4. Instructions. Each instruction occupies one address. Such an instruc- 
tion always contains an operation code for the left-hand six binary digits. 
These two octal digits designate which of the 42 instructions the computer 
is to perform. The right-hand four octal digits usually indicate the address 
from which data are to be taken for use with the instruction being performed. 
Some instructions do not need an associated address; they use the allotted 
digits for other purposes. In the shift instruction, for example, these digits 
indicate the number of digit-positions to be shifted. In any case the third 
and fourth octal digits from the left are not interpreted by the computer 
in the execution of an instruction. 

While the computer is in operation, the arithmetic section operates 
asynchronously with the memory section, each waiting in turn for the other 
to perform its function. In normal operation, when the address of the first 
instruction has been sent to the memory unit (to be held in the Drum 
Address Register, DAR), further action is suspended until that address has 
passed under the playback heads and its contents have been transferred to 
the X register. While the 6-binary-digit operation code is being decoded, 
the 12-binary-digit address, if contained in the instruction, is transferred 
back to the Drum Address Register. No arithmetic operation can be per- 
formed until the specified address has been located on the drum and its 
contents have been transferred to the X register. Actually an address is 
found electronically by circumferentially counting timing channel pulses 
and matching the count with the contents of the DAR. While the instruc- 
tion is being executed, action in the memory unit is suspended. Then the 
address from which the next order is to be taken is transferred to the DAR. 
In this way the computer normally proceeds in linear sequence from one 
instruction to the next. The Instruction Address Register, IAR, holds the 
address of the next instruction in linear sequence to be sent to the DAR 
before the execution of each instruction. In normal sequence ‘‘one”’ is auto- 
matically added to the contents of the IAR as each instruction is carried out. 
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Four “jump” instructions provide important exceptions to the linear 
sequence described. Each one of these can specify the address from which 
the computer is to take its next instruction under a certain condition. 
That is, the four right-hand octal digits of these instructions themselves 
are transferred to the DAR under certain conditions; this transfer interrupts 
the normal transfer of the contents of the IAR to the DAR. These four 
octal digits are also transferred to the IAR, so that operation again proceeds 
in a linear sequence starting with the instruction at this new address. This 
ability to alter the linear sequence of instructions, together with the ability 
to alter instructions by performing simple arithmetic with them, is the basis 
for the real flexibility of this general purpose computer. 

Following is a brief description of the instructions: 


1) Fourteen insert instructions: These cause the transfer of 24 binary 
digits replacing the previous contents of a receiving location. Included are 
transfers between the drum, Accumulator and Q register. Besides these 
normal transfers, provision is made for the transfer of absolute values, 
complements and specified digits, and for the storage of the right-hand 15 
digits (address part) from A. 

2) Eleven arithmetic operations: They include special additions and 
subtractions to facilitate double precision (48-binary-digit) arithmetic and 
manipulations with absolute values, and summation of. the contents of 
successive storage locations. In division the non-negative remainder is re- 
tained in A. Both clear multiply and accumulate multiply produce 48- 
binary-digit products in the Accumulator. 

3) Four jump instructions: One is unconditional. Two others depend 
on the sign of A or Q (control is transferred if the register contains a negative 
number). The fourth is a zero test for the Accumulator. 

4) Two shifting operations: These instructions provide for the circular 
left shifting of the contents of A or Q. 

5) Three stops: optional, intermediate (which depend on pre-set 
switches), and final. 

6) Three output instructions: One causes the typewriter to print a single 
character or perform a single function (carriage return, space, etc.). Another 
causes the typewriter to perform as above and also causes the punching 
of the corresponding teletype code on tape. The third causes the contents 
of successive storage locations to be punched in a form suitable for input 
or later printing. 

7) Five logical operations: These include two digit-by-digit multiplica- 
tion instructions, noncarry addition, clear half of A and a special ‘transfer 
and add one” instruction. 


Division and multiplication (depending on the number of one’s in the 
multiplier) require approximately 2} seconds. All other operations which 
refer to storage can be performed at the rate of 220 per minute, i.e., one 
every two drum revolutions or .272 seconds per single address operation. 
Those which do not require reference to the memory (shifts, jumps, trans- 
fers between A and Q) are executed with twice the speed. An exception to 
these rules is the zero-jump test of the Accumulator. This instruction, 
though reference to the memory is not necessary, requires the time of 2 
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drum revolutions since ‘‘one”’ is automatically subtracted from the contents 
of A after the zero test. 

5. Input, Output, and Operation. Normal perforated input tapes con- 
tain seven channels across the tape, six of which are used for binary infor- 
mation. Four lines, therefore, contain 24 binary digits of information—a 
complete number or instruction. Each line can be considered to have one 
octal digit on each side of the tape feed hole. The seventh channel normally 
contains a punch in every fourth line of the tape to signal the end of a 
24-digit group. For normal input the operator inserts the tape in the reader 
and indicates the address into which the first group is to be stored by pressing 
digit keys of the Drum Address Register. When reading has started, suc- 
cessive groups on tape are automatically stored at successive addresses as 
long as holes continue to appear in the seventh channel in every fourth line. 
This normal tape input takes about 24 second for each 24 binary digits, 
i.e., 90 addresses are loaded per minute. Another type of drum loading from 
tape is possible. In this case each line (6 binary digits) of tape is loaded into 
the least significant six places of consecutive storage locations. This type of 
drum loading proceeds at 220 storage locations per minute. 

Under manual control, not only can the operator inspect in lights the 
contents of any desired memory location, but he can cause insertion of 24 
binary digits there. The contents of A, Q and X are always displayed in 
lights even during computation. Under manual control of the operator, as 
well as when program-controlled, the machine is capable of automatically 
punching the contents of successive storage locations onto tape in a form 
directly suitable for input when the first address is specified by the operator. 
Output of this kind is at the rate of 110 addresses per minute. Also, when 
operator- or program-controlled, the machine can sum the contents of con- 
secutive storage locations. This operation, used for checking the storage of 
a complete program or a set of data whose sum is known, proceeds at the 
rate of 440 addresses per minute, leaving the cumulative sum in the Accu- 
mulator. 

Actually the operator, when performing the manipulations discussed 
above, and the computer itself, while under control of the program, make 
use of three additional 24-binary-digit registers. The Storage Insertion 
Register contains a number to be stored in the memory unit; the Storage 
Output Register holds the number when it comes from the drum; and the 
Storage Blocking Register holds binary ‘‘1’s” in digit positions where stor- 
age of a number is to be blocked. All three may be loaded manually by 
pushing buttons for binary 1’s. 

After both data and instructions are in the memory, the calculation may 
be started (there is no program instruction by which the computer can call 
for additional data to be stored on the drum during computation). When 
starting computation, the operator by means of a switch may cause the 
machine to perform a single instruction or to execute instructions suc- 
cessively from one to the next in normal operation. The computer stops 
automatically and indicates the trouble if the operation code it receives is 
a nonpermissible one, if it reaches a meaningless print code, or if in division 
the quotient is too large for the Q register. 

6. Maintenance and Use. Ideally, after a 15 to 30 minute warm-up 
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period and a successful run of the test program each morning, the computer 
would be available for operation during the remainder of the day, since 
no routine maintenance periods are scheduled. The warm-up period may 
be prolonged or the operation interrupted for needed nonscheduled main- 
tenance. Since no marginal checking features, computational checks or 
indications as to the faulty section are provided internally, much time is 
sometimes consumed in finding and remedying trouble. Ordinarily programs 
are written to include computational checks or are written so that a com- 
plete duplication of each computation section is performed before proceeding 
to the next. 

The Logistics Research Project usually operates the computer 8} hours 
a day during a five day week. For the first eleven months after delivery in 
May ’51, the computer was maintained by several of The George Washing- 
ton University undergraduate student engineers on duty a total of about 
70 hours per week. This indicates that computer time was devoted to edu- 
cation of new maintenance men. Since then, one full-time engineer and one 
half-time undergraduate student have been maintaining the computer. 
During the first year of operation at the University, about $400 was spent 
for replacement parts and approximately $100 for parts for computer 
modification. This includes the cost of replacement of 21 electronic tubes 
and 50 relays. The maintenance men have undertaken projects other than 
routine maintenance. As an example, a test rack has been assembled which 
allows checking of any of the computer’s 48 removable electronic chassis 
apart from the computer. During this first year of operation, approximately 
42.5% of the time during which the machine was on was spent in mainte- 
nance, original checkout and education of maintenance men, and 2% for 
improvement or modification of the computer itself. Thus 55.5% of the 
time was available to the operator-programmers for the year. This last 
figure includes time spent in input, output, program checking and actual 
computation (which sometimes turns into time spent pinning down the 
cause of a machine error). For the first six months of operation the operator- 
programmer time was 43.4% of the total time; this increased to 68.3% for 
the other six months of the year. The computer, then, is purposely a research 
and educational tool in maintenance as well as operation. 

Several problems from organizations outside the Logistics Research 
Project have been allotted machine time. These problems have been dis- 
cussed in advance to determine eligibility for use of the machine with the 
Logistics Branch of the Office of Naval Research, Department of the Navy, 
or with the staff of the Logistics Research Project.' Tape preparation, dupli- 
cation and print-out accessories are available to those whose problems have 
been accepted. 

Following is a selected list of programs which were run during the first 
year, together with pertinent comments: 


1) The Project used the machine to compute quarterly manpower re- 
quirements implied by a proposed four-year shipbuilding schedule. Given 
the building dates and total manpower requirement of a ship, a fourth- 
degree polynomial was evaluated to distribute the man-hour requirements 
into yearly quarters. These requirements were summed for the fleet. 

2) A program was prepared for the iterative solution of games by 
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Brown’s method—each player having 5 strategies. The method provides w 
convergence to the value of the game within .01 after about 8 hours of 
computing (1000 iterations) for the games computed. 


3) Solution of simultaneous linear equations by Crout’s method was T 
performed for m = 5, and a general program was prepared for m < 10. 5 

4) The Project used the computer to assist in obtaining a research w 
solution to a fleet logistics planning problem. 

5) Currently, a calculation is being performed which obtains informa- if 
tion pertaining to feasible schedules for tanker delivery of fuel. Later this p 
information will be used to obtain an optimal schedule with the hope of th 
application to monthly delivery of petroleum products by the Military Sea e 
Transport Service. ce 

J. Jay Woxr 0 


The George Washington University ; 
Washington, D. C. ( 

This work was sponsored by the Office of Naval Research. 

1 Publications available from the Logistics Research Project are: Jos. B. KRUSKAL, 
Jr., A Programmer's Description of ABEL, A Magnetic Drum Relay Computer, Orders for the a 
ONR Relay Computer, a detailed description of the instructions, and ONR Relay Computer 
Operators Manual. (; 


On the Rounding Off of Difference Tables ( 
for Linear Interpolation ( 


In order to simplify linear interpolation many tables contain the differ- 
ences A(x) = f(x + d) — f(x) between two consecutive values of the func- 
tion. Since the values of f(x) given in the tables are usually rounded off ( 
the question arises whether the value of A(x) given in the table must be the 
difference between the rounded off values f(x) of f(x) and f(x + d) of 
f(x + d) or A, the result of the rounding off applied to A(x). It appears on the 
first view plausible that we obtain better results in the second case since we s 
use here more information about f(x). However, the detailed analysis shows 
that this is not so. If the values f(x) are given with m decimals so that the 
rounding off errors do not exceed h = } 10-*, the part of the interpolation 
error due to the rounding off of f(x) and f(x + d) does not exceed h, if the ( 
difference used is f(x + d) — f(x), while if we use A, this error can come 
arbitrarily near to 2h. 

Since this situation is not apparently realized by all computers of tables, ( 
I should like to develop an observation on this subject which was published 
elsewhere.! 

We give first an example. 7 

In computing the decimal logarithm log 9684.8 we start from the values k 
log 9684 = 3.986 054 78; log 9685 = 3.986 099 63 and from the rounded 
off values log 9684 = 3.986 05; log 9685 = 3.986 10 with an error 
<h = 34 10~. 

We have 





A = 4.485-10-, A = 4-10-. 


We obtain then by the ‘‘complete”’ interpolation 
log 9684.8 = 3.986 054 78 + 0.8-4.485-10-> = 3.986 090 66, 
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while in using the rounded off values of both logarithms and A we obtain 
3.986 05 + 0.8-4-10-> = 3.986 082. 


The difference is .866-10-* > h, while, if we take instead of A the difference 
5-10-* of the rounded off values, we obtain the value 3.986 09 and the 
whole rounding off error is .066-10-° < h. 

From the following theoretical discussion it follows in particular, that 
if the rounding off errors of f(x) and f(x + d) can be considered as inde- 
pendent random variables, the probability that the rounding off error of 
the value of f(x + td), obtained in using A, is >(1 + )h, OS <1, is 
equal to 3(¢ — 2n + n?/t), if 1 > ¢ > ». Without loss of generality, we can 
consider a function f(x) with the two values f(0) = fo, f(1) = fi. Then for 
0 <¢# <1, linear interpolation gives 


(1) Lf =(1-Ofot+th = ft td, A=fi-— fo. 


For an h = } 10-* we denote generally by 4 the result of rounding off of 
a to n decimals, so that we have generally 


(2) la—a| =Ah, & = 0 (mod 2A). 
We use the notations 
(3) fi=fite, lel sh, @=0,1), 
(4) A=fi— fo. 
With A we form 
(5) L* = fo + “A, 
(6) L— L* =@+#(A4 — A). 


We consider now for a constant 9 with 0 = n < 1 in the &, «-plane the 
set S(t, 7) of points for which we have 


\L — L*| > (1+ »)h, 


that is, 

(7) & + #4 — A)| > (1+), 
where 

8 A — A = «& — & (mod 2h), 
(8) |A—A| <h. 


Our problem is to determine the area of S(t, 9). The inequality (7) can 
be satisfied only for t > 9, since |@| = h, |A — A] SA. 
Ki We have 
(9) | €o| = k, |e| = aA. 


Further we must have 
le. — &| > h, 


since otherwise A would be f; — fo and in that case the left side of (7) is - 
We will have therefore two cases accordingly as 


lA 
> 


a a4—-a >h or II. 64 —@ < —h. 








In both cases we have 
(10) €1€9 < 0. 


CaseI. hZ=a>O0>gZz=—-ha-—ea>h. 
Here we have from (8) 


A-—-A =a — &© — 2h, 
and (7) becomes 
(11) €) + te — @) — 2th< — (1+ )h, 


since the left hand expression in (11) cannot be >A. 
The points satisfying (11) lie beneath the straight line 


(12) (1—fAo +tat (1+ — 2h =0. 
The line (12) meets « = — h/ in the point 
P(e = —h, a = (1 — 0/t)h) 


and the line «, — ¢ = h in the point 


Pi(€ = (t a oe 1)h, q = (¢ = n)h). 


The corresponding part of the set S(t, 7) is therefore the shaded triangle 


in the diagram 


%, 


h 


1 
(t-%|)h 





-h 


and has the area 3(¢ — »)*h?/t. In the case II we have an entirely symmetric 
argument and have therefore for the total area of S(t, 7) the expression 
(t — »)*h?/t, while all points (¢€, €:) satisfying (9) cover the area 4h. 

We see that the “geometric” probability that the point (é, €,) will fall 
into the set S(t, n) is }(¢ — »)?/t, an expression that is not only positive for 


t > n, but even fairly large. 


A. M. Ostrowski 


American University 
Washington, D. C. 


This work was performed on a National Bureau of Standards contract with the Ameri- 


can University. 


1A. OstrowskKI, Vorlesungen tiber Differential und Integralrechnung. Basel, 1951, v. 2, 


p. 294-296. 
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Mechanical Quadrature Near a Singularity 


The purpose of this note is to present coefficients to facilitate computa- 


b 
tion of integrals of the type f xi f(x) dx. If the integrand x—!f(x) has no 
0 


singularity at the origin, the coefficients presented are still applicable, 
though in this case other well-known procedures for numerical integration 
might be more appropriate. In this regard, it should be noted that suitable 
modifications of f(x) permit one to employ the coefficients in evaluating 


6 
f xf(x)dx. In applied problems, the analytical form of f(x) is usually 
0 


known though it may be rather complicated in form. The singularity may 
be removed by imposing the transformation x = y’ or by fractional inte- 
gration after which any one of a number of numerical procedures can be 
applied. The former approach is not satisfactory since the upper limit of 
integration as a rule becomes irrational and so complicates tabulation of the 
resulting integrand. Integration by parts requires the computation of f’(x). 
If f(x) is quite complicated, its derivative is usually even more so, and again 
an excessive amount of tabulation is required. A third possibility is to ex- 
pand f(x) into a power series about the origin. However, in many cases 
this is not easy owing to the complexity of f(x). Furthermore, the series 
may be slowly convergent. A more satisfactory procedure is to write the 
Lagrangian polynomial which precisely fits f(x) at equally spaced points 
over the interval (0,5). Then a numerical integration formula is easily 


b 
evolved by calculating the moments f x-ixmdx, m = 0,1,2,---. It is 
0 


convenient to consider the integral in the form A, = f xi f(x)dx. Let 
0 


f, = f(rh),r = 0,1, ---, and write 
(1) A, = 2(nk)' (2 fn) D-— + R., 
r=0 


where R, is the remainder term. Following the procedure described above, 
the exact coefficients y,“ and D, are given in Table I for m = 1(1)10. 

It is clear that the (m + 1) point formula so derived is exact if f(x) is a 
polynomial of the mth degree. If f(x) is not of this form, it is of interest to 
examine the remainder. Evaluation of the remainder term in linear methods 
of approximation has been discussed by MILNE.'? Following his notation, 
the remainder is given by 


(2) a f ™ for0(s)G(s)ds. 


The analytical form of G(s) is easily deduced from the approximating 
formula. For the ( + 1) point formula 


(3) n'G(s) = i x—i(x — s)"dx — 2(nh)! > (rh — s)"y-™/Dna, 
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where 


(rh — s)” = (rh — s)"if rh>s 
(4) = 0 ifrh<-s. 


The formula for A, is exact if f(x) = (x — s)" and so 


nh n 
(5) 0 = f x(x — s)"dx — 2(nh)' > (rh — s)"y-™/Dy. 
0 


r=0 
Suppose 0 < s < hk. Subtracting (5) from (3) and using (4), a straight- 
forward reduction gives 
Qntl/ bee 1)"2!s"+3 
1-3-5----- (2n + 1) 





(6) m!G(s) = 2(—1)"(nh)}yo™s" — fO<s<h. 


Denote the right hand side of (5) by Gi(s). Let 


HAs) = 2(nh)*(rh — s)"y-@/Dna; 


(7) G,(s) = Gya(s) + H,-1(s), y = 2,3, -+*, &. 


Then 
(8) n'G(s) = G,(s) if (7 —1h<s<rh, r =1,2,---,n. 


We extend some terminology by Milne and say that G(s) is definite 
if it does not change sign in the interval of (0, 2h); otherwise, it is indefinite. 
If G(s) is definite, application of the mean value theorem to (2) gives 


nh 
(9) R, = porno) f G(s) ds, 


where 0 < 6 < mh. In any event, bounds for the error are given by 


nh 
| < uf |G(s)|ds_ or R,,| < MKnh, 
0 
M = Max | f®*(s)|, K = Max |G(s)|. 


O<s< nh O< s<nh 


(10) 


If G(s) is indefinite, it must vanish at least once in the open interval 
(0, xh) (observe that G(0) = G(nh) = 0). Suppose é is the one and only one 
such point. Then the remainder may be composed of two parts. 


E nh 
(11) R, = porno) f G(s)ds + pron) f G(s)ds 
0 g 


where 0 < 6; < — < 0. < mh. The extension of this argument for G(s) 
vanishing at more than one point in the open interval of integration is 
obvious. 


nh 


It is known that f G(s)ds may be evaluated directly from the approxi- 
0 


mating formula without a knowledge of G(s). In the present instance 


nh n 
(12) f G(s)ds = 2(nh)*h=+{n"+1/(2n + 3) — SS rt y,™/D,}/(n + 1)! 
0 


r=1 


Employing the remainder term for the polynomial approximation to f(x), 
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an expression equivalent to (12) is 
nh nh 
(13) f G(s)ds = f x(x — h)(x — 2h) --- (x — nh)dx/(n + 1)! 
0 0 


Thus, if G(s) is definite, computation of a rigorous error term is considerably 
simplified. As to the definiteness of G(s), no general theorem is available. 
But employing (8), we can compute G(s) for each formula and ascertain 
in a heuristic fashion if it is definite. In this study, G(s) is definite for all 
the even point formulas. For the odd point formulas, G(s) vanishes at one 
point in the open interval (0, zh). For each integration formula given in 


nh 
Table I, the exact value of f G(s)ds is presented in Table II. For the 
0 


£ nh 
odd point formulas, the values of &, f G(s)ds and f G(s)ds are also 
0 g 


tabulated; the first mostly to 5D, the latter to 6D. 
The coefficients y- were checked by verifying that (1) is exact for 


nh 
f(x) = x", m = 0,1, ---,. The values of f G(s)ds were obtained em- 
0 


ploying (12) and checked by (13) and integration of (8). The function G(s) 


TABLE I 
VALUES OF y-™ AND D, 
In each column headed by 2, the first coefficient is ys, the second is y;", etc. The 


n 
last number is the value of D, = 2 y,-™. 
r=0 


n=1 n=2 n= 3 n=4 n= 5§ n=6 

(2-point) (3-point) (4-point) (5-point) (6-point) (7-point) 
2 6 34 250 972 15498 
1 8 45 416 1685 31032 
— 1 18 24 40 —7965 
3 — 8 224 840 26480 

15 — 31 460 — 3870 
105 _— 161 12312 
945 —- 1588 
4158 a 
75075 

n=7 n=8 n=9 n= 10 
(8-point) (9-point) (10-point) (11-point) 

7 66808 59 61306 548 91535 24888 70076 
16 03182 138 56896 1328 43888 65511 43600 
—5 35080 —82 58912 —923 11164 —62119 84725 
15 04055 203 11680 2220 74370 1 55860 46400 
—2 88120 —134 03240 —1681 73334 —1 74234 77400 

5 53602 142 49344 1659 86415 1 97949 48768 

3 48488 —32 57376 —407 45628 —1 26966 12600 

1 01115 44 69632 342 44694 84722 73600 
——— 5 30095 245 60415 — 20655 00900 
40 54050 —_———_ 61 17959 18057 96400 

344 59425 ————_—— 1976 69471 
3394 89150 


1 64991 72690 
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TABLE II 
ERROR COEFFICIENTS 
13 
. G(s)ds/h"\(nh)"2 clk J, G(s)ds/heo J. ™ G(s)ds/hn*32 
1 —2/15 
2 8/315 1.24092 0.039468 —0.003551 
3 —1/35 
4 16/1485 2.69858 0.022193 — 0.000644 
5 —1018/81081 
6 152/25025 4.11461 0.015041 —0.000163 
7 — 35098 /49 22775 
8 180 95776/45831 03525 5.53190 0.011218 —0.000050 
9 —26 20473/5658 15250 
10 1061 15816/3 79480 97187 6.9556 0.008860 —0.000017 


was expanded in a Taylor series about an approximate value of ~, and & 
was determined by inversion of this series. Taylor’s formula was also used 


g 
to calculate f G(s)ds. The latter calculations were performed in duplicate 
0 


and are correct to within one unit of the last decimal place given. 
It appears that coefficients of the kind presented here were first de- 
rived in a thesis by M. Bartes.* This work gives exact coefficients for 


nh 
f x™!2f(x)dx, m = +1, nm = 2(1)4. The remainder terms are examined 
0 


in some detail. In a thesis by M. E. YOUNGBERG,‘ similar formulas were 
given for m = 3(1)7. With the exception of the seven- and eight-point 
formulas (the eight-point formula is totally in error), exact coefficients are 
also tabulated to enable integration over any sub-interval defined by the 
points at which f(x) is tabulated. No estimates of the error are given. 

In a recent paper, E. L. KApLAN® has derived three- and five-point 


rah 
coefficients for the evaluation of f x™!2F(x)dx, m = + 1. All coefficients 
rh 


are in decimal form. For three points, coefficients are given for r; = 0(1)18. 
In addition, some coefficients are presented so that, using three ordinates, 
the integral over the first two or last two ordinates can be found. The five- 
point case is similarly treated, but the values permit integration over 
three ordinates at most. Two sets of coefficients are tabulated. One is exact 
for consecutive powers of x as in our case. The other is exact for f(x) an 
even function. Finally, error terms are given for selected cases having the 
former property. In effect, it is assumed that the mean value theorem is 
applicable and the error coefficients are derived by evaluating equations of 
the type as of the right hand side of (12). But the author rightly remarks 
that the results must be used with caution if rz — r; > 1. As a further aid 
in evaluating the remainder, Kaplan also gives values which account for 
the second approximated term in the Taylor series expansion of f(x). We 
have examined the G(s) function for m = 2, 7, = 4, re = 6 and m = — 1 
and find that it too vanishes in the open interval of integration. The same 
is true for the three-point formula where 7; = 0, re = 2 and m = 1. 
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On the Numerical Solution of Equations 
Involving Differential Operators with 
Constant Coefficients 


1. The General Linear Differential Operator. Consider the differential 
equation of order n 


(1) Ly +> Fy, x) = 0, 
where the operator L is defined by 


Ly = ¥ Pale) k? 


and the functions P;(x) and F(y, x) are such that a solution y and its first 
n derivatives exist in 0 < x < X. In the special case when (1) is linear the 
solution can be completely determined by the well known method of varia- 
tion of parameters when m independent solutions of the associated homo- 
geneous equations are known. Thus for the case when F(y, x) is independent 
of y, the solution of the non-homogeneous equation can be obtained by mere 
quadratures, rather than by laborious stepwise integrations. It does not 
seem to have been observed, however, that even when F(y, x) involves the 
dependent variable y, the numerical integrations can be so arranged that 
the contributions to the integral from the upper limit at each step of the 
integration, at the time when y is still unknown at the upper limit, drop out. 
Thus again the computation can be made to involve merely quadratures. 

It is not often that the solution of the homogeneous equation can be 
simply determined, and it is perhaps for this reason that attention has not 
been given heretofore to the possibility of simplifying the numerical evalua- 
tion of the solution by making use of the solutions to the homogeneous 
equation. However, in the case when the functions P;(x) in L are constants, 
the solution of the homogeneous equation is easy to determine. This is 
particularly true when the order of the differential equation is fairly low. 
In the instance when the operator L is of second order, with constant coeffi- 
cients, the method of using the integral equation often has decided advan- 
tages over the usual methods employed for solving differential equations. 
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For this reason attention will now be centered on a second order operator 
with constant coefficients. 
2. Linear Differential Operators of Second Order with Constant Coeffi- 
cients. Let Ly now be specialized as follows: 
d*y dy 
(2) Ly = be +b , + cy, 


where b,c are real constants. Let y(0) = po, and y’(0) = p; be assigned. 
The differential equation (1) can be replaced by the integral equation 


(3) y = aye™™ + ave™ — G(x) + G(x), 


where 
a + a2 = po; am, + ame = pi, 


Gi.(x) = Jom f e~mkt Fy, na|/ — 4c)!; k =1,2, 
0 


and m, mz are the roots of m? + bm +c = 0, provided b? — 4c ¥ 0. In 
the special case when b? — 4c = 0, (3) becomes 


(3a) y= oa + ax + f (t — x)e*F(y, nat], 


% =p, @ = pt po. 


It should be observed that when b? — 4c is negative, m, and mz, are 
conjugate complex numbers. When 3, c, and the initial values are real, the 
imaginary component of (3) will drop out, and the discussion which is to 
follow will apply to this case and to (3a) as well. For the sake of simplicity, 
therefore, we shall now assume that m, and mz are real and distinct. Let 


Xe =X t+rh, x =X +sh, t =X + ph, yy, = y(x-). 
With the above, G;(x) takes the form 


(4) Gi(x) = em2-20G,(x9) + C(x), 

where 

(5) Ci(x) = em *h(b? — 40) f e—mkPh Fy, x9 + ph)dp. 
0 


If the integrand of (5) is approximated by an (s + 1)-point polynomial, 
then 


(6) f e~mkrh Fd pb _ > a,e—™rh Fy, Xr) + Ri, 
0 


r=0 


where the coefficients a, result from the integration of the polynomial and 
the truncation term, R;, can be represented by 


(7) R, -f $(t)[t, Xo, x1, +++, Xe dt. 
0 


In (7) [t, x0, x1, --+, Xs] is the divided difference! of order (s + 1) of the 
function e~-™ F(y, t) and ¢(t) is the polynomial approximation of e~™F(y, ¢). 
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Thus 
(8) — Gi(x) + Go(x) = — e™@-G\(x9) + € 2 Go(x) 


s--1 
+ h(b? — 4c)+ ¥ aF(y,, x, fermm™e-? — em} — Ri + Reo. 
r=0 


Note that the term involving F(y,, x.) dropped out, since the coefficient of 
this term appeared with the same sign in G, and G2. Although this fact is 
well known in the theory of integral equations, its importance from the 
viewpoint of the numerical evaluation of the solution needs emphasis. Thus 
the evaluation of — G; + G2 does not depend on the value of the function at 
the end of the interval. We therefore do not need a “predictor’”’ formula 
(using MILNE’s terminology”). The steps of integration at any stage can 
therefore be carried out as follows: 


1) Evaluate — G;(x) + G(x) by (8). 

2) Compute a; exp (mx) + a2 exp (mex); hence knowing (— G; + Gs), 
we now know y,, from (3). 

3) Knowing y,, we can use (6) to evaluate the integral G,(x) by a mere 
quadrature. G2(x) is now also known, since G;(x) and — G,(x) + Ge(x) 
are known. 


3. The Truncation Term. If m, and m, are negative, e~™'F may require 
a higher order approximation formula than F(¢) itself (although by no means 
necessarily so). In any case, # must be small enough so that R, in (7) is 
indeed negligible for the accuracy aimed at. In some cases it is actually 
possible to take a /arger step h when the integral equation is used, than the 
one that can be taken when the differential form (2) is operated with. 
Moreover, the process of solving the integral equation may be more stable 
than the corresponding solution of the differential equation by stepwise 
integration. This is especially true when (5? — 4c)! is large numerically. 
The example following illustrates the case. 

4. Example. Consider the differential equation 


(9) u’ +vuu' +e*+ f(u) =0, 


where 


B 
f(u) = exp (4 _ 5) : 


This differential equation occurred in connection with certain steady state 
solutions, and among the various parameters which were used, one set was 
the following: 


A = 74.997736; B = 257.42325; d = 2.19885; (0) = 1.294; 
u'(0) = — M = — 10; vo = 24.38. 


Suppose we attempted the evaluation of (9) by the usual method of 
stepwise integration, first computing an integral of w’’, using (9), then 
integrating u’ to obtain u. An analysis of the manner in which an error in 
u’ is propagated shows that the interval h would have to be taken no larger 
than 1/(5v) over the entire range of the integration. Thus for v = 25, an 
interval as small as 0.008 would be needed. 
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A study of the behavior of f(u), which depends not only on v but also 
on M and u(0), shows that for very small values of x, it is sometimes neces- 
sary to take # even smaller than 1/(5v). This is true whether we solve the 
equation by integrating (9) or by using (3). However, when M = 10, for 
example, and the integral equation is used, it is possible to use a step h as 
large as 0.1, for x larger than 0.4, and yet we would have to maintain an 
interval of about 0.008 if the differential equation were used. It can be 
verified that it is possible to lose all significance in u’ within a relatively 
short range of x when (9) is used, unless the interval h is kept small enough. 
When (3) is used, on the other hand, the solution remains very stable, and 
the size of the interval can be increased just as quickly as the behavior of 
f(u) permits. It turns out that e*f(u~) behaves no worse than f(x) itself. 
The successive derivatives of f(u) with respect to x are numerically very 
large near the origin, but they go down in magnitude to reasonable levels 
at x = 0.4, where M = 10. For this problem, therefore, the form (3) is far 
superior to the usual method of solving the differential equation. 

The permissible magnitude of / near the origin, when (3) is employed, 
can be determined by usual methods, and the necessary starting values can 
be computed in a simple manner; hence it is not deemed worth while to give 
the specific details of the solution. We shall, however, re-write (3) to apply 
specifically to (9). Here 


F = f(u) + e*; m =0; m. =—v; po = u(0); p = —M. 


When the above are substituted into (3), the equation becomes 


1 = 1 x 
u=Ui-~ f flu(a)at + ee [ enptu( iat, 











where 
wa Ur M 1 
U; = - +5 | a - ]+~-4* , wv— 1) #9. 
v—1 v v—1 v 

Thus let 

i: £e 1 t ; ge 

G(x) -<f fdt; S(x) = wee f ef dt — ~{ f dt, 
) Jo v 0 v Jo 


H(x) -1e f "ef dt = S(x) + G(x). 


v 


Using the notation 


?: = $(sh) 


for any function ¢, we have, by the previous analysis, 


1 1 
Ui = Ue41 + - eS. —1 + G,-1 | — 0 Ge + Ps+t) 


1 (v+1)h 1 (v+1)h 
don = ervey f ef dt — - f f dt. 
( 


v—1)h US (w—1)h 


where 


If @.41 is evaluated by Simpson’s rule, we have 


h 
bers = 35 [HE — 1) fa + (CM — Ifa] + R+ Suan 
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where 6,4; is the rounding error. The total error at x can be approximated, 
roughly, by 


1 Ist of 1 :. of 
———_— = peal. — g-rsti vt — 
€s41 -f (6+ R) ar dt + Je ae (6+ R) a dt. 
Differences can be used in the usual manner to approximate the magnitude 
of derivatives. The recommendations for the special example are: 


a) Compute ¢,4: by Simpson’s rule or some other suitable integration 
formula. 

b) Evaluate 1,41. 

c) Evaluate G,,; by quadrature (since u,4; is now known to the required 
accuracy). 

d) Write Sya1 = Usaa — Uses. 


Return to (a) for the evaluation of ¢,42 in the succeeding step. 

The process can be readily coded for high-speed machines. 

One may inquire whether (3) is always the better form to use, compared 
with (1), or whether there are ranges of the parameters where one of the 
forms is better than the other. An examination of the way in which the 
constants enter into the solution shows that when b? — 4c is large form (3) 
would always have advantages over (1), since the square root of this quan- 
tity enters into the denominator of some of the terms. In the case when 
this quantity is small, however, it is likely that (1) would be the better form 
to operate with. 

5. The General Case. It may be worth while to remark that the method 
applies to the general case mentioned in Section I, where the coefficients 
P,(x) are functions of x. The solution can be constructed formally by using 
the usual method of variation of parameters for linear equations.* The 
complete solution will be of the form 


n 
y=> ui, 
k=1 


where the m functions u; are the known solutions of the homogeneous equa- 
tion, and the m integrals V, must be generated numerically by using the 
Wronskian of the solutions u;. In the numerical process, the value of y at 
some point x is obtained from the property that the contribution to y(x) 
from the upper limit of integration drops out. After y(x) is known, each 
individual V;(x), which is required for carrying forward the numerical 
process, is then evaluated by a mere quadrature. 


GERTRUDE BLANCH 
National Bureau of Standards 


Los Angeles 24 
The preparation of this paper was sponsored (in part) by the Office of Naval Research, 
SN. 
1 See L. M. MiLne-TuHomson, The Calculus of Finite Differences, London, 1933. 


2W. E. MILNE, Numerical Calculus, Princeton, 1949. 
3E. L. Ince, Ordinary Differential Equations, London, 1927, section 5.23. 
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RECENT MATHEMATICAL TABLES 


1018[D].—J. Peters, Sechsstellige Werte der Kreis- und Evolventen-Funk- 
ttonen von Hundertstel zu Hundertstel des Grades nebst einigen Hilfstafeln 
fiir die Zahnradtechnik. Zweite verbesserte Auflage, Bonn, Ferd. Dumm- 
lers Verlag, 1951, viii, 222 p., 14 X 20.8 cm. 


The first edition of this work (viii, 182 p.) appeared in 1937 and the size 
of its page was 19.1 X 26 cm. For the second edition a photographic copy 
of these pages was made in reduced size, the only changes being in the 
addition on page viii of a preface for the second edition, and of nine new 
words on the title page. It will be observed, however, that an Appendix 
of 40 new pages has been added. These pages are mainly filled with tables 
of results useful in studying problems dealing with teeth of gears. 

In the main body of the volume, with a° as argument at interval .01, 
one may read off on opposite pages 6D values of the 6 trigonometric func- 
tions; of the arc a; of the polar coordinates 6° and sec a of the evolute of a 
unit circle; of the radius of curvature, tan a of the evolute; and of arc 4 
= eva = tana — arca. 


R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


1019[F, L].—L. A. DRAGONETTE, “Some asymptotic formulae for the mock 
theta series of Ramanujan,’”’ Amer. Math. Soc., Trans., v. 72, 1952, 
p. 474-500. 


The author is concerned chiefly with the function 
fM=UVMA+*it?@---a+tgq}- 
n=0 


introduced by RAMANUJAN. This function generates the function A(m) 
defined by 


f@ = y A(n)q", 


which is similar to the partition function p(m) in that A(m) possesses an 
asymptotic expansion similar to that of Harpy & RAMANUJAN.' Table 1 
(p. 495) gives the exact values of A(m) for » = 0(1)100. The asymptotic 
formula is 


A(n) = 2 dal) (h(n — 1/24))-t exp {(4/k)(1/6)(m — 1/24)*} + O(n! log n). 
k<n 


The coefficients \,(2) are real and given in Table II (p. 497) as trigonometric 
polynomials in for k = 1(1)14. 
D. H. L. 


1G. H. Harpy & S. Ramanujan, “Asymptotic formulae in combinatory analysis,” 
London Math. Soc., Proc., s. 2, v. 17, 1918, p. 75-115. 


1020[G, K].—F. N. Davin & M. G. KENDALL, ‘‘Tables of symmetric func- 
tions. Parts II and III,” Biometrika, v. 38, 1951, p. 435-462. 


Part I of this paper appeared in 1949 and is reviewed in MTAC, v 4, 
p. 146, where also will be found references to other symmetric function 
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tables. The tables of Part I gave expressions for the monomial symmetric 
functions 
(a1, a2, ---) = 2) xi%x2™ --- 


in terms of the power sums 
a = (r) _ Xv xy" 


and inverse tables for all weights <12. 
Part II gives the monomial symmetric functions in terms of the elemen- 
tary symmetric function 
a, = (1°) = > xx -- + x, 
with inverse tables. 
Part III gives the monomial symmetric functions in terms of the homo- 
geneous product sums h, defined by the generator . 


1+ bit + hf + --- = (1 — at + af — ---), 


with inverse tables. 

Parts II and III extend also to weights <12 and have the same format 
as Part I. The tables of Part II have been compared with earlier tables. 
The tables of Part III are new. 

D. H. L. 


1021(G].—M. Osma, “On the irreducible representations of the symmetric 

group,” Canadian Jn. Math., v. 4, 1952, p. 381-384. 

A table is given of the number of irreducible representations of the 
symmetric and alternating groups of order m for m = 2(1)40. If a represen- 
tation by a Young diagram has its rows and columns interchanged we obtain 
an associated representation. The number of self-associated representations 


is also given. 
Dt 


1022[I].—Z. Korat, P. Carrus, & K. E. KAvANaAGu, “A new formula for 
repeated mechanical quadratures,” Jn. Math. Phys., v. 30, 1951, p. 
44-48. 


The authors start with the Hermite interpolation formula 


(1) fe) = h(x) f(a;) + E fila) f(a) + [pa(x) PFO” (u)/(2n)! 


where 

hj(x) = {1 — (x — as) [pn (as)/Pn' (as) HUA) F, 

hj(x) = (x — a,)[I,(x) }, 

I(x) = pa(x)/[(x — @5)Pn'(a;)), 

Pa(x) = (x — a)(x — ag) --+ (x — ay), 
and w is an inner point of the range including all a; and x. If the coefficients 
a; are determined so as to satisfy 


1 
f hjx)dx =0, j =1,2,-+++,m, 
a1 


then the well-known Gauss quadrature formula results. The authors de- 
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rived new formulas by setting 

1 1 
(2) f hy(x)dx = f hj(x)dx =0, j =2,3,---,m: 

me —l 
Let 

1 
60) =F); fi fede = 2F/(a) + U, 


1 z 
U = fief o(y)dy. 


When equations (2) are satisfied, there results the quadrature formula 


where 


(3) U = 2f(a:) + H,o(a;)) +R 


where the remainder, R, involves the factor ¢°"-")(y); hence R = 0 if $(y) 
is a polynomial of degree no higher than 2” — 2. The authors tabulate the 
weight factors H;“” to 6, 7, or 8D, and the points a; to 7D for n = 2, 3, 4, 
and to 6D for m = 5 and 7. They discovered the interesting fact that the 
system (2) yields more than one set of real solutions a; for the values of n 
considered in sufficient detail. The following examples for » = 3 illustrate 
the character of the coefficients: 


Type A (a, = 0) Type B (a; ¥ 0) 
a, =0 a, = +.5889711 
@2= .5477225 Hz =  .3042903 dz = +.2250452 He = F.8445196 
@3 = —.5477225 H; = —.3042903 a3 = +.5293628 Hs; = #.3334225 


No attempt was made by the authors to prove the existence of solutions of 
(2) for general »; but for odd values of n, solutions of Type A can be shown 
to exist. 

The new formulas have advantages and disadvantages similar to those 
of the Gauss quadrature formula. It is of course true, as the authors infer, 
that by a simple linear transformation U can be made to represent the more 


general double integral 
d x 
W -f ax f o(y)dy. 
e b 


1 y 
W=U,4+V; Ui =3(d- of aw f oi(u)du; 
-1 a 


Indeed 


Vi = Hd — 0)? f” alway, 
where 
hed = 0 +0) +26 ~ ah: 2 hc - OG - 0. 


However, it must be remembered that both V; and U; need to be evaluated, 
and that the argument for which ¢ must be computed may be cumbersome 
in practice. Where the calculations are not prohibitive, a further simplifi- 
cation is in fact possible. For U vanishes whenever ¢(u) is an even function 
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of wu. Hence if ¢(—y) is defined, write in (2) 
o(y) = 3(¢3 + bs); os = O(y) + O(—¥); os = O(y) — O(—9). 


It is necessary to apply the quadrature formula only to ¢,, if a formula of 
Type A is used, and thus only half the number of multiplications have to 
be performed. In many cases arising in practice ¢, may be just as simple to 
evaluate as ¢ itself. 

There is a misprint in the authors’ formula (21); in the coefficient of the 
cosine term replace the term 2x? of the numerator by 2°. 


GERTRUDE BLANCH 
NBSINA 


1023[K].—H. A. Davin, ‘‘Further applications of range to the analysis of 
variance,’’ Biometrika, v. 38, 1951, pp. 393-409. 


Consider the mean of & uncorrelated ranges for samples of size » from 
normal populations with common variance o’. Let ¢ be the expected value 
of [(mean range)/o ?. Let » equal the number of degrees of freedom for the 
x’-distribution which is approximately equivalent to the distribution of 
(mean range)/co. Table I contains values of c and » for m = 2(1)10 and 
k = 1(1)5, 10. Table II contains values of c and » for use in a double classi- 
fication analysis. Here m = 2(1)9 and k = 2(1)10, 20. In Table III a split- 
plot design with / main treatments, m blocks, and N subtreatments (N is 
reviewer’s notation) is considered. Table III contains values of ¢ and 
vy’ = v/l for m = 2(1)10 and N = 2(1)9. Values of ¢ are given to 2D, those 
of vy and y’ to 1D. Let d, be the expected value of (range)/o for a sample of 
size n while V, is the variance of this quantity. Table IV contains values 
of dn, dn/ Vn, dn?/Vn to 2D for use in analyses of single classification with 
unequal cell frequencies. Here m = 2(1)20. 

J. E. WaALsH 
U. S. Naval Ordnance Test Station 
China Lake, California 


1024[K].—P. M. Grunpy, ‘“‘The expected frequencies in a sample of an 
animal population in which the abundances of species are log-normally 
distributed. Part I,”” Biometrika, v. 38, 1951, p. 427-434. 


Given N species in a population with the expected number of individuals 
per species (abundance), m, following a log normal distribution 


f(m)dm = 





e~ (in (m/a)]?/20% dy, | 


1 
V2xmo 


where o” is the variance and In a the mean of In m (hence a is the median 
abundance). The probability of obtaining r members of a given species in 
a sample is e~"m"/r!. Hence the expected proportion of species having r in 
the sample, ¢,, is 





_ f ” £m pln) den. 


r! 
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Table 1 presents values of 1 — ¢o (expected proportion of species in the 
sample) and Table 2 presents values of ¢: (expected proportion of species 
with singletons) to 4D for o? = 2(1)16 and logisa = — 2.0(.25)3.0. 


R. L. ANDERSON 
North Carolina State College 


Raleigh, North Carolina 


1025[K].—-P. G. Guest, “The fitting of polynomials by the method of 
weighted grouping,” Annals Math. Stat., v. 22, 1951, p. 537-548. 


Given equally spaced values of x and the corresponding 2 observed 
p 

values of y, it is desired to fit a polynomial u,(x) = } b,x! (pb < n) by a 
j=0 


method of weighted grouping. To obtain an unbiased estimate of by, = dp, 
the points are divided into not more than 2p + 1 groups of successive 
points, the number of points in the ith and in the (m — 7 + 1)th groups 
being equal. The sum of the y’s in each of p + 1 of these groups is to be 
assigned a non-zero weight. Because of symmetry, not more than (p + 2)/2 
different weights are involved if is even. If pis odd, p + 1 different weights 
are involved but half of these are the negatives of the other half. The 
groupings are determined for each m and # in such a way that the variance 
of a», assuming equal variances of the y’s, is a minimum. To estimate ,,;, 
j =0,1,---,p—1, the a;, 7 = 0,1,---,p—1 are calculated and the 
estimates of b,; are obtained from the relation, estimate of b,; = @; + Bj41, j@j41 
+ --+ + Bay. (Alternate f’s are zero.) 

A table (p. 541-545) gives the groupings, the weights (exact), the f’s 
(Bo and Bs, exact; Be and B53, 10S; Bo and Bs, 9S) for m = 7(1)55, 
j = 0(1)5. 

The relative (compared with least square procedure) efficiencies of the 
estimates of b,; are discussed and a table of some limiting relative efficiencies 
is given. The lowest limiting relative efficiency appearing in the table is 
.889 for estimates of a;. Computational schedules and an example are given. 

K. J. ARNOLD 
University of Wisconsin 
Madison, Wisconsin 


1026[K].—H. O. Hartley & E. S. PEarson, ‘‘Moment constants for the 


distribution of range in normal samples,” Biometrika, v. 38, 1951, p. 
463-464. 


The moments indicated by the title are the mean y,;’; the central 
moments pr, 7 = 2(1)6; o = V2; Bi = ms*/pe*; Bo = wa/ms?; Ke = wa — 3yn?; 
Ks = os — 10usue; Ke = we — 15p4u2 — 10p3? + 30u2*. These moments are 
tabulated for m = 2(1)20 except for us, us, x4, Ks, Ks, Which are given for 
m = 2(1)12. The tabled values are given to the number of decimals followed 
in parentheses by the number of significant digits as follows: u;’, 5(6); 
Mea, 5(S); 3, 4(4); Ma, 3(4) ; Ms, 3(3); He, 1(2-3); Co, 4(4); Bi, 4(4); Bo, 3(4) ; 
4, 3(2—3); xs, 2(0—2); xe, 2(0-3). 

W. J. Dixon 
University of Oregon 
Eugene, Oregon 
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1027[K].—L. KAARSEMAKER & A. VAN WIJNGAARDEN, Tables for Use in 
Rank Correlation, Report R 73, Computation Department of the Mathe- 
matical Centre, Amsterdam, 1952, 17 p. mimeographed. 21.0 X 32.4 cm. 


For a given number m of subjects ranked according to two criteria, the 
size of KENDALL’s coefficient of rank correlation is proportional to the 
score S.' If S is the sample value, on the null hypothesis of no association 
E(S) = 0; to test this null hypothesis P,(S = S), the probability that in 
n subjects S = S, on the assumption of no association, is used. It is assumed 
for this paper that no ties in ranks occur. Kendall? tabulated P,(.S = S) 
for m = 1(1)10 for all possible values of S and the authors extend this table 
giving P,(S =S) to 3D for m = 1(1)40. They divide the table into two 
parts, Tables I and II, for S even and odd, respectively. Table III gives for 
n = 4(1)40 and for a = .005, .01, .025, .05 and .1, the smallest value of S 
for which P,(S 2S) < a. Since for > 40 the normal distribution is a 
good approximation to that of S, values of the standard deviation of S, os 
are listed in Table IV for m = 40(1)100 to 5D. There is a typographical 
error in the formula for og where it is first written on p. 2; it appears cor- 
rectly on p. 3 and in the table. A recurrence relation used in the calculation 
of Tables I and II is developed. 

.. ©. Se 

1M. G. KENDALL, Rank Correlation Methods, London, Griffin, 1948. 

2 Loc. cit., p. 141. 


1028[K].—K. C. S. Prtxat, “Some notes on ordered samples from a normal 
population,”’ Sankhyd, v. 11, 1951, p. 23-28. 


Consider a sample of size from a normal population with zero mean. 
Table I contains 1% and 5% points of T = 2 (midrange)/(range) to 2 and 
3D for m = 3(1)10. Table II presents values of coefficients for approximate 
determination of the distribution of the median (average of two central 
order statistics) to 6D for m = 4,6,8. Next consider two independent 
samples from normal populations. Form the range for each sample and let 
the sample size be m, for the larger range and m2 for the smaller range. 
Table III contains 1% and 5% points of F’ = (larger range)/(smaller 
range) to 2D for m, m_. = 2(1)8. Let o; be the standard deviation for the 
population yielding the larger range and a2 the standard deviation of the 
other population. Table IV contains a power function comparison of the 
F’-test and the Snedecor F-test to 3D for (m, m2) = (3,7), (4,4), (8, 2), 
(5, 2), (8, 3), and o:/o2 = 1(.5)2.5. 

J. E. WALsH 
U. S. Naval Ordnance Test Station 
China Lake, California 


1029[K].—J. E. Watsu, ‘Some nonparametric tests of whether the largest 
observations of a set are too large or too small,’”’ Annals Math. Stat., 
v. 21, 1950, p. 583-592. 


This paper is apparentiy the first published work covering nonpara- 
metric tests for the detection of outlying observations. For the nonpara- 
metric tests, the usual restriction of normality can, of course, be removed 
and the sample considered to be drawn randomly from one or more con- 
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tinuous symmetrical populations. Assuming that the sample values are 
arranged in increasing order of magnitude: x(1) < x(2) < --- < x(m), the 
tests considered are (a) the detection of whether the r largest observations 
are too large to be consistent with the hypothesis that the populations 
from which the sample values came have a common median, (b) whether 
the r largest observations are too small and (c) whether the populations 
are symmetrical in the tails. With regard to (a) and (b), the proposed tests 
also cover similar hypotheses for the r smallest observations because of 
symmetry conditions. The tests, based on the order statistics, covered in 
the paper have the somewhat fascinating property in that the Type I error 
or significance level, a, for the sample statistics used is independent of the 
sample size for values of m permitted. (Generally speaking, m should be 
large, in which case the significance level tends toward the value a; how- 
ever, for no admissible value of does the significance level exceed 2a.) 
The details of the actual tests are rather complicated for description 
here, but Table 1 gives the necessary specifications for 21 tests and the 
corresponding significance levels, a, to 4D, for use with r > 4. 
F. E. GRUBBS 


Ballistic Research Laboratories 
Aberdeen Proving Ground, Md. 


1030[L].—S. CHANDRASEKHAR, DONNA ELBERT, & ANN FRANKLIN, “The 
X- and Y-functions for isotropic scattering. I,” Astrophys. Jn., v. 115, 
1952, p. 244-268. 

S. CHANDRASEKHAR & DONNA ELBERT, ‘‘The X- and Y-functions for 
isotropic scattering. II,’’ Astrophys. Jn., v. 115, 1952, p. 269-278. 


The functions in question are solutions of the integral equation 


X(u) = 1+ poo f (w+ LX (w)X(u’) — V(u) Y(u’) Jay’, 


V(u) = e7 + pus f (u — CY W)Xw’) — X(u) Vw) Ty’ 


Approximations were described in CHANDRASEKHAR’S book Radiative Trans- 
fer, Chapter VIII. 

In the first of these papers the approximate solutions were improved by 
iteration and are tabulated for 0 < » < 1 for the following values of the 
parameters: wo = .95, .9, .8, .5, 7 = .05(.05).25, .5, 1. The tables are to 4 
or 5D, and the difference between the tabulated solution and the corrected 
second approximation (of Chandrasekhar’s book) is also given. 

RIn the second paper, ‘‘standard solutions” are defined in the case w) = 1. 
They satisfy the relations 


[xa =2, [rua =o. 


The moments of X, Y, X,, Y; are called an, Bn, an*, Bn*. The laws of diffuse 
reflection_are expressed in terms of the solutions 

X*(u) = X.(u) + QulX.(u) + Y.(u)], 

Y*(u) = Y.(u) — QulX.(u) + Y.(u)], 
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where 
Q wien a; — 6° 
(ay + Bi)r + 2(a2 + B2) : 
Table 1a gives ao, Bo, ao*, a, 81, and 
§ = 1 — {(2 — woao)ay + woBoh:} 
for 7 = .05(.05).25, .5, 1 and w) = .95, .90, .80, .5. Table 1b gives ay + Bo, 
a + Bi, a2 + Bo, ar", Bi, (ay*)? — (B;*}?, O, and 3 for the same values of + 
and wo = 1. 
Table 2 gives X*, Y*, X,, Y, for wo = 1, 7 = .05(.05).25, .5, 1, 
u = 0(.01)1. 
Table 3 contains certain combinations of the moments. 





A. E. 


1031[L].—W. CueEsTER, “The reflection of a transient pulse by a parabolic 
cylinder and a paraboloid of revolution,” Quart. Jn. Mech. Appl. Math., 
v. 5, 1952, p. 196-205. 


Table I (p. 204) gives 4D values of 


x2 4 r 2 —1 
ar f gt iE ~ re f eas] + re*| dr 
0 0 


oo @ 2 =f 
+f e-(rtie {[: + $re7 riens| + tere*| dr 
0 


—T 


and of 


for r = 0, .1, .2(.2)2. In the last expression the integral with respect to z is 
meant as a Cauchy principal value. 
A. E. 


1032(L].—G. Diemer & H. DijKGraar, ‘‘Langmuir’s £, 7 tables for the 
exponential region of the J, — V, characteristic,”’ Philips Res. Rep., 
v. 7, 1952, p. 45-53. 


” 
The integral f [f(x) + m*}-idx, where 


f(x) = e7(1 + erf x!) — 1 — 2(x/x)!, 
was tabulated by FREEMAN! in the range 0 < 9 < 20,0 < m? < 20. Modern 
microwave diodes and triodes make it desirable to extend these tables 
considerably. The present paper contains 3S values for 10-? < » < 60, 
0 < m* < 10°. The intervals vary. The two tables overlap and there are 
discrepancies in the overlapping regions. 
A. E. 


1J. J. FREEMAN, “Noise spectrum of a diode with a retarding field,” NBS, Jn. of 
Research, v. 42, 1949, p. 75-88. 


1033[L].—S. GotpstTErn, ‘“‘On diffusion by discontinuous movements, and 
on the telegraph equation,” Quart. Jn. Mech. Appl. Math., v. 4, 1951, 
p. 129-156. 


The y(m, v) may be defined by the recurrence relation (r + 1)y(" + 1, ») 
= r[y(n, » — 1) + y(n,» + 1)] — (r — 1)y(m — 1, »), together with the 
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initial conditions 
y(i,1) = y(1, —1) =4, y(1,7) =Oif»y ¥ +1. 


The author obtains the generating function of the y(n, v), also integral rep- 
resentations and an asymptotic expansion for m large and v? = O(n). To 
test the asymptotic expansion, he tabulates (p. 137) 7D values of (15, v) 
for r = 3,1, 2,5 and » = 1(2)15. AE 


“ 1034(L].—HaArRVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 
23: Tables of the error function and of its first twenty derivatives. Cam- 
bridge, Mass., Harvard University Press, 1952, xxviii, 276 p., 19.5 X 
26.7 cm. $8.00. 


The functions tabulated in this volume are 


1 ’ d” 1 
(-—1) = “nee —t? /2 (n) an “it an oe ~ 
ci) (x) | e dt, ? (x) =| ated » n ama 


The ¢™ are connected with Hermite polynomials. 


Tables of these functions are listed in section 15 of the FMR Index. 
Additional tables are mentioned in MTAC, v. 3, p. 521, and tables on 
punched cards are listed in MTAC, v. 5, p. 203, 204. 

Modern developments in statistics, probability theory, mathematical 
physics, and noise and communications theory made new tables desirable 
over a more extensive range with a finer interval, and, most of all, including 
derivatives of higher orders. The authors of the present work expect that 
their tables will satisfy all present needs. 

Tables | and II give 6D values for ¢—” to @, where x runs from 0 
at intervals of .004 up to the point where ¢ ceases to change. The limiting 
value is .5 for 6, and zero for all other ¢“. The last x increases with , 
being 4.892 for nm = — 1, and 8.236 for m = 10. 

Tables III and IV give values to 7S or 6D for ¢@" to ¢, and x runs 
from 0 at intervals of .002 up to the point where ¢ vanishes to 6D iden- 
tically. This point is 8.518 for m = 11, and 10.902 for n = 20. 

The introduction contains a chapter (by WARREN L. SEMON) giving a 
collection of formulas and integral representations for ¢‘") and their connec- 
tion with Hermite polynomials; a brief section (by the same author) 
explaining the computation of the tables and the steps taken to eliminate 
errors; a chapter (by DAvip MIDDLETON) on applications of the functions 
tabulated in this volume; and a 10D table of all the zeros of all the functions 
tabulated in this volume. Each section is accompanied by a list of references. 

The tables maintain the high standards set by previous volumes of this 
series. AE 


1035(L].—Tu. LarsLe, ‘‘Héhenkarte des Fehlerintegrals,’’ Zeit. angew. 
Math. Physik, v. 2, 1951, p. 484-486. 
Relief diagrams of erf (x + iy) for (i) O< x < 5,0 < y <6, and (ii) 
0<x<1.7, 0 < y < 2.2. 4S values of the real and imaginary parts of 
the first 5 zeros (in the first quadrant) of erf (z). AE 
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1036[L].—C.-B. Linc, ‘Tables of values of the integrals f x™dx/sinh?x 
0 


and f x™dx/cosh?x,”” Jn. Math. Phys., v. 31, 1952, p. 58-62. 
0 


This paper gives values of the above integrals multiplied by the auxiliary 
factor p™*!/2?(m)! mostly to 5D for all admissible pairs of » and m, 
bp = 1(1)8, m = 0(1)15. The reciprocal of the above factor is given exactly 
or to 6S for the same range. C. W. NEtson' has tabulated the first integral 
to 12D for m = p = 1(1)40. Entries common to both tables are in perfect 
agreement. 


YUDELL L. LUKE 
Midwest Research Institute 


Kansas City, Missouri 


1 C. W. Netson, “A Fourier integral solution for the plane-stress problem of a circular 
ring with concentrated radial loads,” Jn. Appl. Mech., v. 18, 1951, p. 173-182. 


1037[L].—R. C. Lock, ‘‘The velocity distribution in the laminar boundary 
layer between parallel streams,”’ Quart. Jn. Mech. Appl. Math., v. 4, 
1951, p. 42-63. 


The “‘standard solution,’ g(¢), of the Blasius equation ff’ + 2f’”’ =0 
is that solution which has the expansion 


g(é) = —1+ eH — def + pee 
for — < 0. 
Table II (p. 58) gives 4D values of g, g’, g” for = — «, —16(2) — 10(1) 
— 5(.4).2, and .5537 = &, where g(&) = 0. 
A second solution, g:(¢), defined for — > &, depends on a parameter pu. 
It is a solution of the Blasius equation, subject to the initial conditions 


gi(f) = 0, gi’(f) = g’(&), a1'"(fo) = (pu)*g’’ (£0). 
With a = [g,'(~)]}-!, b = &, two further functions f; and f2 are defined by 
film) = agi(am + 6), fe(n2) = ag(ane + 5). 


These are also solutions of the Blasius equation, m > 0, nz < 0. 

Tables III-VI (p. 59-62) give 4D values of fi, fi’, fi’, fo, fo’, fe” for, 
respectively, pu = 5.965 X 104, 100, 10, 1. 7 and m2 range over their respec- 
tive ranges at intervals suited to the physical problem. 

Table VII (p. 63) gives, for the case py = 1, a similar tabulation of that 
solution for which f2’(— ©) = .5014 (instead of 0 as in the previous cases). 


A. E. 


1038[L].—B. Misura, ‘‘Wave functions for excited states of mercury and 
potassium,’’ Cambridge Phil. Soc., Proc., v. 48, 1952, p. 511-515. 


_ The radial wave function P(n/|r) for excited state (, 1) of Hg satisfies 
the differential equation 


d? 2Z,(r) i+ 1) 
a + (nl) + - 7 





P(nl\r) = 0, 
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where r is the radial distance, ¢(n/l) is an energy parameter, and Z,(r) is the 
total effective nuclear charge. The wave function is normalized so that 


in [P(nl|r)}*dr = 1. 


For the numerical integration the variables p = In (10007) and S = r~iP 
were used. 

Table 1, p. 512, gives 3 to 5S values of Z, and 4D values of S for the 
states (6s), (6p), (6d), (7s), (7p), (7d), for p = 0(1/3)2(1/6)4(1/12)11-25. 
The values of Z, were derived from existing tables for the neutral atom, 
and the wave functions found by numerical integration of the differential 
equation. 

The differential equation for K contains an additional term V,(r) 
= — ta(r? + r?)-*. Table 3 gives 4D values of P(4p|r) for r = 0(.02) 
.3(.05).6(.1)1.2(.2)3.2(.4)16(1)30. 

Tables 2 and 4 give the numerical values of certain integrals. 

A. E. 


1039[L].—A. Papou.is, “On the accumulation of errors in the numerical 
solution of differential equations,” Jn. Appl. Phys., v. 23, 1952, p. 
173-176. 


By solving a differential equation with various initial conditions the 
author exhibits a method for determining the error in the numerical solution 
of such equations. Consider the differential equation (f) y’ = f(x,y) and 
assume that the interval [a, b] is divided into m subintervals by the points 
of subdivision x;. Let 6, be the difference between the value of the correct 
solution at x = x, and the computed solution y,; at x;. (8; is the sum of 
truncation and round-off errors.) If we assume that y(x) is correct up to 
Xx-1 and from x; to b satisfies (t), then at x = b we shall have an error 
Ex = y(b) — yn. Let y(x) be a solution of (ft) with the initial condition 
(x0) = yo and Y(x) a solution of (¢) with V(x) = yo + M (M aconstant); 
then under the assumption (A): r(x) = Y(x) — y(x) is small compared 
with y(x), the author shows that the total error R is the sum of the E; 
( = BxLr(b)/r(xx) ]), and hence for x large, 


R = [mr(b)/(b — aif [B(x)/r(x) Jdx. 


Under assumption (A), f(x, Y(x)) — f(x, y(x)) = r(x)e(x), where g(x) 
= Of(x, y(x))/dx and hence r(x) is determined from the differential equa- 
tion r’(x) = o(x)r(x). If B(x) is the result of round-off errors alone and if 
they are random with variance d, then the dispersion of R is 


[a?nr*(0)/(6 — a)} f [1/r%¢e) Jae. 


Extensions to systems of two equations, dx/dt = f(x, y, t), dy/dt = g(x, y,t), 
are given in detail with some remarks on systems of higher order. 


K. S. MILLER 
New York University 


New York, N. Y. 
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1040[L].—A. VAN WIJNGAARDEN, Table of the integral 


1 
f exp (—v~? — xv)v~?dv. 
0 
Mathematical Centre, Amsterdam, Report R 176, 1952, 6 mimeographed 
leaves. 


The integral mentioned in the title is given to 5D (“the last decimal is 
not absolutely reliable’’) for 


x = 0(.25)3.5(.5)14.5, 
p = 0(1)10. 


For x = 0 the integral becomes 


io) 
1 


an incomplete gamma function. A special table for x = 0 to 6D is given for 


p = — 25(1)11. Several methods were used to compute the tables including 
power series and recurrence relations resulting from integration by parts. 
oo: fe. 


1041[L].—E. M. Witson, Solutions of the Equations (y'’)? = yy' and Two 
Other Equations, Admiralty Research Laboratory, Teddington, Middle- 
sex, November, 1951. 


The first part of this report is concerned with the solution of y” = —Vyy’ 
such that y(0) = 0, y’(0) = 1. Solution depends on the inversion of the 


v 
integral x = f (1 — ##)-dt, where x is the independent variable. A gen- 
0 


eralization of this integral has been studied by R. GRAMMEL (MTAC, v. 5, 
p. 155). When x approaches 4nV3/9, y approaches unity while y’ and y” 
approach zero. Starting at this point the function tabulated is the solution 
of y” = + Vyy’ which admits of an integral representation similar to the 
previous one. Series solutions are also given. A photostat table is available 
giving y to 6D for x = 0(.002)6 and log y to 6D for x = 6(.01)7. The report 
tabulates y to 6D for x = 0(.05).5(.1)6. Author believes maximum error 
is less than 0.7 unit in the last figure given. Second differences, mostly 
modified, are also tabulated. To facilitate interpolation near the origin, 
y = 4x!/15 is tabulated to 6D for x = 0(.05).25. 
In the second part, values of the integral 


8B 
f(B, p) = (2e-*/B*Vm) f To(2pn)e-*ndn 


are tabulated to 4D for 6 = 0(.25)4, p = 0(.25)5. Ip is the modified Bessel 
function. The author claims entries are correct to within one unit of the 
last figure. The table has been subtabulated so that entries are in intervals 
of .05 and is available in photostat form. A more extensive tabulation of a 
function simply related to the above has been made by S. R. BRINKLEY, JR. 
& R. F. BrinkLtey (MTAC, v. 2, p. 221) and S. R. BRINKLEY, Jr., H. E. 
Epwarps & R. W. Smitu, Jr. (MTAC, v. 6, p. 40). 
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Part three tabulates to 3D that zero of usin x — cosx + e~“* which 
lies between + and 2x for u = .1(.01).3(.02)2 and Vu = 0(.02).5. The 
author states that the last figure should be correct to within 0.7 of a unit. 
Linear interpolation yields full accuracy and first differences are provided. 


YuDELL L. LUKE 
Midwest Research Institute 
Kansas City, Missouri 


MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to errata in RMT 1022, 1032. 
214.—E. P. Apams & R. L. HippisLey, Tables of Elliptic Functions, Smith- 
sonian Miscellaneous Collections, v. 74, no. 1, Washington 1939, 1947. 


The heading of page 294 


for E’ = 1.5629622295 
read E’ = 1.5631622295. 


H. J. HAUER 
Box 342 
China Lake, California 


215.—(1) A. A. GERSHUN, ‘‘Berechnung des Volumleuchtens,”’ Physikalische 
Z. d. Sowjetunion, v. 2, 1932, p. 149-185 [MTAC, v. 2, p. 191]. 
(2) E. Scumipt, “Die Berechnung der Strahlung von Gasraumen,” 
Zs. Verein Deutscher Ingenieure, v. 77, 1933, p. 1162-1164. 
(3) S. GotpstEtn, ‘On the vortex theory of screw propellers,’’ Roy. 
Soc., London, Proc., v. 123A, 1929, p. 440-465. 


Three tables are given in (1), on p. 172, 175, and 180, respectively. 
Table I, containing the function 
@ gt 
F(x) = 1 — (1 — xe* + x*Bi(-x) =1- 20 | Sat 

to 4D for x = 0(.01).02, .05, .1, .2(.2)1, 1.6(.4)2.4, 3, 4, was read against 
the same function given in (2) in complementary form on p. 1163, also to 
4D mainly, for x = 0(.01).02, .05, .1(.1)1.0(.2)2.0, 2.4, 2.5, 3(1)5. The 
discrepancies, and the extra values given in (2), were checked, revealing 
the following errors in (1) and (2): 


(1) BF for read 
0.4 4925 4854 
3.0 .9822 .9821 

(2) 0.3 .6000 .6001 
0.9 .2516 .2514 
1.2 .1680 .1679 
1.4 .1296 .1292 
1.6 1011 .0998 
1.8 .0777 .0774 
25 .0328 .0326 
4.0 .00545 .00552 


.00175 .00176 
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Table II, of the function 


© (—1)Ptige 2 = 
Pls) = 2D Ga yy -1- al - A + xe), 
was recomputed. The following values are in error: 
F, for read 
2.8 .8039 .8038 
4.0 8863 8864 
5.0 .9289 .9232 


x/2 
Table III gives F, = 1 — f e-78t cosidt = 4a[Ii(x) — Li(x)], 
0 


where J, and L, are Bessel and Struve functions of imaginary argument. 
Comparison with (3) indicated one error in (1) 


x= for 9281 read 9271 

and two in (3) 
p = 0.8 for 0656 read 0671 
1.2 for 0039 read 0045 


DAPHNE P. KILNER 
Scientific Computing Service, Ltd. 
London, England 


216.—P. R. E. JAHNKE & F. Empe, Tables of Higher Functions. 4th (revised) 
ed. Leipzig, 1948. 


This edition gives on p. 9 a more extensive table of the maxima and 
minima of the sine integral, that is, si (x) for x = #(x)24x. Comparing these 
values with those obtained from interpolating in the NBS tables we find 
the seven errata: 


wy for read 
18 — 0.007673 — 0.017673 
19 +0.006744 +0.016744 
20 — 0.005907 — 0.015907 
21 +0.005151 +0.015151 
22 — 0.004463 — 0.014463 
23 +0.003834 +0.013834 
24 — 0.008258 — 0.013258 
H. J. Gray, Jr. 


Moore School of Electrical Engineering 


University of Pennsylvania 
Philadelphia, Penn. 


On p. 239 the entry immediately to the right of 6.3 in the x column 
should be shifted down one line. Thus —A,“(26.30) = 0.096170. The first 
significant figure of the function does not change from 6 to 5 until x is 
between 6.32 and 6.33. 


J. H. Hancock 





1745 Lanier Pl., N.W. 
Washington 9, D. C. 
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217.—NILs PippinG, “Die Goldbachsche Vermutung und der Goldbach- 
Vinogradowsche Satz,”’ Abo, Finland, Akad., Acta Math. Phys., v. 1i, 
no. 4, 1938, p. 1-25. 


* for m, read 
6944 61 37 
10006 149 83 
23926 47 17 
31004 73 67 


S. A. JOFFE 
515 West 110th Street 


New York 25 


UNPUBLISHED MATHEMATICAL TABLES 
In this issue there is a reference to an unpublished table in RMT 1041. 


150[F].—D. D. WaLL, Table of Wilson’s Quotient. 11 leaves tabulated from 
punched cards. Deposited in the UMT Fite. 


For each of the 709 primes p < 5381 the table gives the least positive 
remainder on division of {(p — 1)! + 1}/p by p. This remainder is zero for 
p = 5, 13, and 563. The table was produced on the IBM Card Programmed 
Calculator. [See also MTAC, v. 5, p. 81, MTE 182. ] 


D. D. WALL 
IBM Corporation 


Los Angeles, California 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
Cannon, 415 South Building, National Bureau of Standards, Washington 25, D. C. 


DISCUSSIONS 
ASYNCHRONOUS SIGNALS IN DIGITAL COMPUTERS 


It is frequently necessary, during the operation of a digital computer, 
to inject signals from sources that are not synchronized with the computer 
itself, for example, the manual signals. This operation may be initiated by 
pressing an appropriate push button. In this discussion, we will not be 
concerned with such problems as ‘“‘bounce’’ of contacts, wavering pressure 
or the possibility of repeated operation because of completion of computa- 
tion before the button is released, but only with the fact that the contact 
is made (or broken) at a random moment with respect to the computer 
timing pulses or “‘clock.’’ Probably the most important source of auto- 
matically generated signals asynchronous with the computer proper is the 
input equipment. Whether data are introduced by magnetic tape, punched 
cards, manual keyboard or other means, it is generally introduced at a 
much lower rate than transfers within the computer itself and at intervals 
which do not synchronize with the main ‘‘clock.”’ 
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The presence of asynchronous signals creates a problem of special type 
in that their effect under certain conditions is not strictly determined in the 
digital sense. That is, examination of such a signal at a time synchronous 
with the machine proper may not be interpretable as either “yes” or “‘no”’ 
over a narrow range of relative timing of the signal to the machine, but only 
as “‘maybe.”’ If an entity similar to one of Maxwell’s famous demons were 
available to convert such “‘maybe’s”’ to either “‘yes’’ or “‘no,”” even in arbi- 
trary fashion, the problem would be solved. Automatic equipment, how- 
ever, responds in continuous fashion to signal amplitudes and durations, 
and these must have a transition range between any discrete set of states. 
When one is dealing with wholly synchronous signals, such transitions can 
be made to take place between examining times and thus do not cause 
difficulty. With asynchronous signals, the random relative timing causes 
randomness in the transition durations and thus permits conditions under 
which examination of the signal may occur during the period of transition. 

One might reply to the above discussion that the existence of an indefi- 
nite ‘“‘maybe’’ is not of much importance since, no matter whether the 
remaining equipment behaves as if it were “‘yes’”’ or as if it were “no,” 
operation would be satisfactory. This is true if al/ the remaining equipment 
made the same interpretation of the ‘‘maybe,”’ and, in fact, in certain cases 
it suffices to ‘“weight’’ one element deliberately so that its interpretation 
of the ‘‘maybe”’ can never be “‘yes’”’ when any of the other elements behave 
as if the “‘maybe”’ is “‘no.’”’ In most cases, however, error results if any two 
elements do not interpret the signal in the same manner. The “maybe” 
condition thus usually implies error. An otherwise perfect machine therefore 
is not error free in the presence of asynchronous signals, but has a finite 
probability of error. The problem is to make this probability extremely small. 

This problem was recognized at an early stage in computer development. 
Thus the designers of the ENIAC incorporated additional “flip-flops” as 
buffers between the asynchronous signals and the remainder of the machine. 
The designers’ reasoning as expressed in lectures and operating manuals 
was essentially as follows. The asynchronous signal, being random, could 
occur at a time when it caused passage of only a partial ‘‘clock’’ pulse thus 
making operation indefinite. By using this partial pulse, however, to set a 
flip-flop, and gating a later pulse through it, a “full” output pulse is ob- 
tained if the flip-flop sets and no output pulse if it remains unset. In the 
latter case, a full pulse will be obtained during the next cycle (since the 
asynchronous signal is of long duration) which will be certain to set the 
flip-flop. This scheme has, apparently, yielded equipment which functions 
quite satisfactorily. There is, however, a flaw in the above reasoning in the 
assumption that the flip-flop is definitely either set or not set when it gates 
the interrogating pulse. From our previous discussion, there must be a 
finite, although very small, probability that the gated pulse is also “partial.” 
It seems obvious, nevertheless, that the method is effective in markedly 
reducing the probability of error. 

This use of ‘trigger’ circuits between an asynchronous signal and the 
remainder of the machine is, to the author’s knowledge, resorted to by all 
computer designers. The purpose of this paper is to point out that the effect 
is to reduce but not eliminate probability of error, and that material im- 
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provement is possible by consideration of the various factors involved in 
such an arrangement. It is easy to see that the region of indefiniteness is 
decreased by “squaring” up of the signal and interrogating pulses, by 
speeding up of the “flipping”’ of the trigger circuit and by increasing the 
time interval between the pulse gated by the signal and acting upon the 
trigger and the pulse gated by the trigger and applied to the balance of the 
machine. The ‘‘squaring”’ of the signal and pulses reduces the probability 
of error in a smaller proportion than the decrease in rise or fall times since 
duration is a factor as well as slope. Moreover, the rise time cannot be 
readily increased indefinitely because of the usual circuit limitations. The 
improvement that can be achieved by careful attention to the shape of the 
pulse is thus quite limited. As concerns time of response of triggers and time 
interval between interrogating pulses, the probability of error will, in most 
cases, be an exponential function of the ratio of the two times. Decrease of 
this ratio by increasing the rate of response of the trigger or by increasing 
the time interval between gating and reading the trigger is therefore ex- 
tremely effective in reducing probability of error. 

We will conclude by indicating, for a particular configuration, how 
probability of error may be evaluated, at least as concerns order of magni- 
tude. For this purpose, we assume perfectly ‘‘square’’ pulses and signal and 
a trigger whose cutput builds up exponentially to a steady state pro- 
portional to the input, and which is provided with positive ‘‘feed-back”’ 
for outputs in excess of both the input and a fixed noise suppression volt- 
age. Such a trigger is approximated by a linear amplifier with ideal diode 
gating. If we take the input as &;, and the trigger output as E, we have 
E = AE,(1i — e~“*), where A is the steady state amplification, ¢ is the 
duration of the input and & is the time constant of the circuit. If E ever 
reaches the value £;, the feed-back replaces the signal and the trigger is 
fully set. If E does not reach the value E) which is the noise suppression 
voltage, no feed-back takes place when the input signal ends and the trigger 
is fully reset. We are interested in the intermediate condition where, at 
time ¢ when the input ends, £; + Ey > E > &). In this case, the feed-back 
becomes effective at time ¢ and we get the equation 


dE 
kat E = A(E - E) 


whose solution is 


Pi jana — emt) — aS. Jewors + Fair 





for the period following cessation of the first interrogating pulse, where T 
is the time interval during this period. Since A > 1 for satisfactory trigger 
operation, this solution represents an exponential build-up. If E falls be- 
tween certain limits, say E, and E;, at the time the trigger is examined, 
error results. The probability of error is thus the probability that ¢ is such 
as to make 
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so that 4; < t < t, where 
th = —kin(S — E./AE,) + (A — 1)T 


and 
te = —kin(S — E;/AE,) + (A — 1)T, 
where 
Eo Ey 
- > eee (A—1) T/k —— 
s-|1 aa. |* +@— DE’ 


If Q is the repetition period for the interrogating pulses, this probability is 
given by 


(tz — th) 
Q 
_k ; S — E,/AE, 
~Q"S—E,/AE,’ 
In view of the natural choice of T large compared to k and E, large compared 
to Eo, S is given approximately by e4-"7* and P is approximated as 
(E; — E,) 
ASE, 


P = 


k 
P=-—Inj1+ 
om 
_ k (Es — Es) 
' Q AE, 
which summarizes the important factors. If, in a particular case, we have 


(E; — E,.)/E, = 0.1; A =2; k =1 microsecond; T = 10 microseconds 
and Q = 25 microseconds, we find 


—(A—1 
e (A TIE 


1 
—— —7 
500 ° , 
which is too high for comfort. However, increase of T or decrease of k by a 
factor of 2 makes 
1 
= — gS = 12 
500 © ‘xi. 
which may be acceptable. Similar evaluation of effectiveness may be made 
for other types of trigger circuits. 


S. LUBKIN 
Electronic Computer Corporation 


265 Butler Street 
Brooklyn 17, N. Y. 
BIBLIOGRAPHY Z-—XXI 


1. Anon., “Analog-digital converter,” Rev. Sci. Instr., v. 22, Oct. 1951, 
p. 800. 


Expository article. 
2. Anon., “Digital computers,” The Electrician, v. 148, March 14, 1952, 
p. 818. 


This is an expository article which very briefly describes certain features 
of the digital computer at Manchester University, England. This machine 
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is one of a series of machines which Ferranti, Ltd., is producing. The second 
one of this series will be located in Toronto, Canada, at the university. 


Epitu T. Norris 
NBSMDL 


3. R. C. M. Barnes, E. H. Cooke-YarsporouGcH & D. G. A. THOMAs, 
“An electronic digital computor using cold cathode counting tubes for 
storage,’ Electronic Engineering, v. 23, Aug. 1951, p. 286-291, and 
Sept. 1951, p. 341-343. 


In this sequence controlled computor, whose numbers contain eight 
decimal digits and a sign, speed was sacrificed to achieve small size and 
simplicity by using relays for switching and ‘‘Dekatron’’ 10-position cold 
cathode discharge tubes for up to 90 words of storage. Numbers and the 
instructions in sequence are read from perforated tape by up to eight tape 
readers and results are printed or punched by up to eight printers or per- 
forators. The computor contains about 380 relays, 18 Dekatrons, 80 therm- 
ionic tubes, and 40 cold cathode triodes, plus 28 relays and 90 Dekatrons 
per 10 words of storage. It occupies three 7 ft. relay racks, together with 
an additional rack per 49 words of storage, one smaller rack for power 
supplies, and a table for tape readers and printers. Total power consumption 
is about 1 kilowatt. 

The computor operates in parallel mode, and for transfer within the 
machine a decimal digit is represented by a train of 0 to 9 pulses in its 
appropriate channel while on the tapes a 2 out of 5 code is used. Round off 
is performed by adding 1 or 0 at random in the 8th place. Negative numbers 
are represented by 9’s complements with end around carry correction. 


R. D. ELBourRN 
NBSCML 


4. H. J. GeIs_er, “R. F. bursts actuate gas tube switch,” Electronics, v. 25, 
Feb. 1952, p. 104-105. 


The article gives a brief description of a gas tube gate, requiring radio- 
frequency excitation of one of its inputs. The r-f energy is coupled to the 
tube through the d-c electrodes and a conducting band around the tube 
envelope. Cold open-circuit resistance is several megohms, and the capaci- 
tance between electrodes and also between the band and the electrodes is 
less than 1 upf. 

The characteristics of the IBM-36 developmental gas tube switch are 
given, and the application of the tube in computer storage and accumulator 
read-in circuits is discussed. 


E. W. C. 
NBSMDL 


5. J. A. Goetz & A. W. Brooke, “Electron tube experience in computing 
equipment,” Electrical Engineering, v. 70, Feb. 1952, p. 154-157. 


The IBM Corporation has 2,500,000 electron tube sockets which are 
used in commercial computation equipment. Two years ago the IBM Tube 
Laboratories established a defective tube analysis program. This paper 
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describes the methods used for preventive maintenance testing of tubes 
before installation, and for analyses of causes of failure or rejection for 
several tube types used in number by IBM. It shows how an appreciable 
gain in machine reliability has been accomplished since this program has 
been in effect. The analysis of common causes of tube failures and vacuum 
tube life expectancy and survival data is of special interest to the computer 
design and maintenance engineer. 


P. D. SHUPE 
NBSCML 


6. W. E. Mutter, “Improved cathode-ray tube for application in Williams’ 
memory system,”’ Electrical Engineering, v. 71, Apr. 1952, p. 352-356. 


This article describes the IBM-79 (RTMA 3VP1), a cathode-ray tube 
designed specifically for use in a Williams’ electrostatic storage system. 
Attempts were made to reduce spot size, to improve “‘spill’”’ characteristics 
and to reduce deflection defocusing and noise pickup. As in all tubes, the 
tube was designed to effect a compromise of conflicting requirements. The 
three-inch size represents a compromise between “bits” per tube and “‘bits”’ 
per unit volume. 

Compared with the 3KP1, the best commercially available tube of the 
same size, there is a decided improvement as shown by the following: 
1) spot separation for a given “‘spill’’ number is about 60 percent of that 
of the 3KP1 in a two-dot test ; 2) spot size is about 76 percent of that of the 
3KP1; and 3) deflection defocusing (from the illustrations) is about 90 per- 
cent of that of the 3KP1 at one inch from the center of the tube. These 
improvements were obtained by reducing the magnification ratio, masking 
down the beam, and shaping the deflection plates. In addition, the capacities 
of the deflection plates, especially to the grid, were balanced to reduce 
deflection caused by beam pulsing. Extreme cleanliness in assembly and 
processing and “sparking” were used to produce a larger percentage of 
blemish-free tubes. Noise pickup may be reduced by an external silver 
coating, grounded near the output end. The smaller spot size gives a re- 
duced, but still ample, output signal. No attempt was made to reduce the 
“gentle rain” of secondary electrons on the storage surface. 


D. C. FRIEDMAN 
NBSCML 


7. T. J. Rey, “On the background of pulse-coded computers, part 1,” 
Electronic Engineering, v. 24, Jan. 1952, p. 28-32; part 2, Electronic 
Engineering, v. 24, Feb. 1952, p. 66-69. 


Part 1 develops the background of the present day digital computer. 
Then the functions of and relations between the Input-Output, Control, 
Memory and Arithmetic units is explained, and the distinction is made 
between serial and parallel operations. The remainder and major portion 
of part 1 is concerned with the representation of numbers within a digital 
machine based on the binary number system and the congruence relation- 
ship of number theory. 

Part 2 describes how the basic building blocks of digital computers may 
be represented by Boolean algebra. There are several examples of the design 
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of such sub-units as an adder from a truth table to a logical circuit. The 
article is concluded by a comparison between digital and analog systems. 


WILuiiaM A. Notz 
NBSCML 


8. Joun J. WILD, “High-speed printer for computers and communications,”’ 
Electronics, v. 25, May 1952, p. 116-120. 


A novel printer, which has produced, under experimental operation, up 
to 900 eighty-character lines a minute without serious degradation of im- 
pression, has been developed by the Potter Instrument Company. The 
device consists of a ‘Flying Typewriter” and electronic control, the former 
being a rotating-type-wheel typewriter and the latter essentially a flexible 
electronic counter with certain special features incorporated to trigger 80 
solenoid-controlled hammers facing the periphery of the type wheel. The 
time when a hammer strikes the type wheel through ribbon and paper 
determines the character printed in the corresponding column. 

The control system automatically distributes the printing of the char- 
acters in the proper order during a revolution of the type wheel, although 
the order of printing is not necessarily in positional sequence around the 
wheel. 

The hammers are resilient to cause them to bounce back swiftly after 
striking the type. Time of contact of hammer and type face is under 0.1 
millisecond, and the time of operation of the hammer is about 2.5 milli- 
seconds. The pulses operating the hammers are timed with a lead of 2.5 
milliseconds to cause the striking of the type character in the center. Even 
at extremely high rates of operation no appreciable blur of print occurs. 
The printing is equivalent, in clearness, to that of a good typewriter. 

The time of operation of each of the 80 hammers is controlled by infor- 
mation set into a trigger tube storage called a PASS unit (printer actuator 
serial storage unit). A motor-driven notched disc and phototube arrange- 
ment provide the necessary synchronizing pulses. The PASS unit has 80 
columns of 6 binary-digit storage, with shifting circuits for loading and 
special gating circuits for driving the printer. PASS may be loaded a column 
at a time, with a 6 binary-digit code, or it may be loaded a pulse at a time. 
Each of its 80 columns can be used as a scale of 10 or 64 counter. In fact, 
the PASS unit is essentially an accumulator and shift register and is capable 
of use by an automatic computer as a part of its arithmetic unit, when 
printing is not being performed. 

The loading of the PASS unit with the information to be printed may 
be thought of as the pre-setting of the 80 counters in the unit. The pulses 
generated by the photo-electric disc are fed individually to step the counters, 
with time lag between columns to prevent the call for printing of the same 
character in more than one column at the same time. When each of the 
various counter columns in the PASS overflows, an output pulse is generated 
which causes the firing of a thyraton and the consequent energizing of the 
solenoid-controlled hammer. The characters around the rotating-type wheel 
and the photo-electric disc are so placed that when the hammer strikes, the 
correct type character is opposite it. 





at 
pul 


cor 
the 
the 


col 


NE 


ha: 
tor 
an 


fol 


ne 


iS. 


= 


—- wt i 


| 


EO a —Ctié«‘SSG 





AUTOMATIC COMPUTING MACHINERY 245 


The new Potter Instrument Company printer was designed to operate 
at high speed, and 2!so to be capable of accepting data for printing from 
punched cards, tapes, or directly from automatic digital computers. Pre- 
setting of selected columns in the control storage inhibits printing in the 
corresponding columns of a line; this feature facilitates the arrangement of 
the output format. This feature, together with the high printing speed, if 
the printer is reliable, should make the unit a useful addition to high-speed 
computer accessories. 

E. W. C. 
NBSMDL 
NEws 


Institute for Advanced Study.—A high-speed electronic digital computing machine 
has been completed and put into operation at the Institute for Advanced Study in Prince- 
ton, New Jersey. The machine is designed to perform very high-speed calculations in pure 
and applied mathematics and in mathematical physics. 

Prior to this public announcement the machine had successfully completed a number of 
quite extensive and important problems. These problems include among others the 
following: 


1) A large number-theoretical problem to test a conjecture, which has never been 
established, of the famous 19th century mathematician, E. E. Kummer. This calcu- 
lation required the instrument to perform about 20,000,000 multiplications and took 
six hours of continuous computing. 

2) Two considerably shorter astro-physical problems, each requiring the solution of 
three simultaneous differential equations. 

3) Solutions of a number of cubic diophantine equations. 

4) Several twelve-hour meteorological predications covering the continental United 
States, each amounting to about 800,000 multiplications and requiring about one hour 
of continuous computing time. This last work is only the first step in an extensive 
research program in theoretical meteorology being carried out at the Institute for 
Advanced Study. 


This machine took six years to design, develop and construct by a staff under the 
direction of Prof. JoHN voN NEUMANN at the Institute for Advanced Study. The initial 
sponsorship of the project came from the Research and Development Service of the Ord- 
nance Corps, U. S. Army. It has continued under the joint sponsorship of that agency, 
together with the Office of Naval Research, U. S. Navy; the U. S. Air Force; and the U. S. 
Atomic Energy Commission. Throughout its history, the project had the support and 
encouragement of the Institute for Advanced Study. In addition, the Office of Naval Re- 
search, since 1946, and since 1951 in cooperation with the Geophysics Research Division, 
Air Force Cambridge Research Center, Cambridge, Mass., has sponsored the work in 
dynamic meteorology. The Office of Naval Research also supported a complementary 
research program in numerical analysis. 

The machine has been the prototype for various subsequent machine developments, 
including three for the AEC and one recently completed by the University of Illinois for 
the U. S. Army Ordnance Corps, Ballistic Research Laboratories at the Aberdeen Proving 
Ground. 

The engineering design was due to JULIAN BIGELOW, its execution to Julian Bigelow 
and JAMES POMERENE, assisted by GERALD Estrin, HEwitt CRANE, RICHARD MELVILLE, 
NorMAN EmMsLIz, and EPHRAIM FREI, as well as by Gorpon KENT, PETER PANAGOs, and 
others. The mathematical and logical design was due to John von Neumann and HERMAN 
H. Go.pstinE. The work on meteorology is under the joint direction of JULE CHARNEY 
and John von Neumann. 
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Data can be introduced into the machine in decimal or in binary form, but the instrument 
proper carries out the calculation in the binary number system since the use of this number 
system is more convenient electronically. Each number handled by the machine consists 
of a sign and 39 binary digits which is the equivalent of a decimal number with a sign and 
approximately 12 decimal places. The machine produces 2,000 multiplications per second, 
1,200 divisions per second or 100,000 additions per second. For a machine of this degree of 
precision it is the fastest one now operating. The machine consists of four principal organs: 
namely, an arithmetic organ which carries out the processes indicated above, a control 
organ which executes the instructions given the machine, a memory organ in which both 
the numerical data of the problem and the instructions which characterize the problem 
are stored and, lastly, an input-output organ which intervenes between the human operator 
and the machine itself. The memory is a system of 40 cathode ray tubes, based on an 
invention of F. C. WiLLtaMs in Manchester, England. The machine can get access to the 
memory organ in 25 microseconds. This organ is capable of storing 1,024 numbers each of 
40 binary digits. 

It is unusually small in physical size, being approximately 2 X 8 X 8 feet, and its total 
power requirements including ventilators, etc., are about 15 kilowatts. It contains about 
2,340 vacuum tubes, almost all of which are double triodes. 

Reeves Instrument Corporation.—Project Cyclone Symposium II on Simulation and 
Computing Techniques was held in New York City, April 29-May 2, 1952, under the 
sponsorship of the Reeves Instrument Corporation with the approval of the U. S. Navy 
Special Devices Center. The program consisted of three sessions, under the chairmanship 
of RAwLEy D. McCoy, Reeves Instrument Corporation. 


Session I 
Analogue Computer Techniques and Appli- | Wednesday, April 30, Morning: 
cations 
REAC solution of problems in structural C. W. BRENNER, Mass. Inst. of Technology 
dynamics 
The use of an analogue computer and G. MartTIN and L. M. LeGatskI, University 
feedback theory for the solution of of Michigan 


structural problems in the static case 
Application of the Electronic Differential G.M.Corcos, R. M. Howe, L. L. Raucn 
Analyzer to eigenvalue problems and J. R. Seviars, University of Michi- 
gan 
Some REAC techniques employed at the L. Pope, David Taylor Model Basin 
David Taylor Model Basin 
Afternoon: 
Simulation studies of a relay servomech- N. P. Tomitnson, Goodyear Aircraft Corp. 
anism 
Use of the REAC asa curve fitting device C. H. Murpuy, Ballistic Research Labora- 
tories, Aberdeen Proving Ground 
Precision in high-speed electronic differ- H. BELL, JR., and V. C. Ripeout, Univer- 


ential analyzers sity of Wisconsin 
On an application of the use of analogue J. H.LANING, JR., and R. H. Battin, Mas- 
computers to methods of statistical sachusetts Inst. of Technology 
analysis 
Solution of linear differential equations C.E. Howe, R. M. Howe and L. L. Raucs, 
with time varying coefficients by the University of Michigan 
Electronic. Differential Analyzer 
Session II 
General Papers Thursday, May 1, Morning: 


JAINCOMP computers and their appli- D.H. Jacoss, The Jacobs Instrument Co. 
cation to simulation problems 
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The Decimal Digital Differential Ana- 
lyzer CRC 105 as a tool for simulation 
and checking analogue computer solu- 
tions 

Problems encountered in the operation of 
the MIT Flight Simulator 

Automatic REAC operation for statistical 
studies 


Mathematical error analysis for continu- 
ous computers 
Checking analogue computer solutions 


Analogue computation of blade designs 


Solution of partial integral-differential 
equations of electron dynamics using 
analogue computers with storage de- 
vices 

Session III 
Computer Components 
Modification of REAC 


REAC servo response 

Applications of differential relays to solu- 
tion of REAC problems 

The design and test of a Linear Wiener 
Filter 


The role of diodes in an Electronic Differ- 
ential Analyzer 
A high accuracy time division multiplier 


An AM-FM Electronic Analogue Multi- 
plier 
Discussion Period on Electronic Multipli- 
cation 


E. Wetss, Computer Research Corp. 


W. W. SEIFERT and H. Jacoss, Jr., Massa- 
chusetts Institute of Technology 

R. R. Bennett and A. S. Futton, Hughes 
Aircraft Co. 

Afternoon: 

F. J. Murray, Columbia University 


W. F. Ricumonp, JR., and B. D. Loveman, 
The Glenn L. Martin Co. 

D. B. Breepon, M. M. MatrTHews and 
E. L. Harper, Westinghouse Electric 
Corporation 

C. C. WanG, Sperry Gyroscope Co. 


Friday, May 2, Morning: 

J. W. Foxutn, Jr., G. F. Emca and F. M. 
Watters, Applied Physics Laboratory, 
Johns Hopkins University 

A. H. MILier, University of Minnesota 

L. M. WarsHAwskKyY and W. Braun, Wright 
Air Development Center 

G. Nestor, Avion Instrument Corp. 


Afternoon: 

C. D. Morritt and R. V. Baum, Goodyear 
Aircraft Corp. 

E. A. GoLpBEeRrG, Radio Corporation of 
America 

W. A. McCoo , Naval Research Laboratory 


UNIVAC Acceptance Tests.—On April 22, 1952, the third UNIVAC system, this one 





constructed for the U. S. Army Map Service by the Eckert-Mauchly division of Remington- 
Rand, Inc., under NBS contract, passed the final test for its acceptance. This was the 
same magnetic tape reading and writing test given the second UNIVAC (see MTAC, v. 6, 
Apr. 1952, p. 119). It required reading 142 million decimal digits and writing 85 million; 
however, in this test the central computer controlled 10 tape systems. Error detecting 
circuits stopped the computer in 8 out of 25 fifteen-minute test units; however, there were 
no undetected errors, and the operator in every case corrected the trouble without assistance 
from an engineer or maintenance technician. That this rate of stoppages, barely low enough 
for acceptance, was not achieved until the seventh attempt at the test indicates that a 
considerable improvement in reliability of tape reading would be desirable. However, this 
should not be used for unfavorable comparison with any other machine. To this writer’s 
knowledge no other machine approaches the UNIVAC’s speed in reading and writing on 
magnetic tapes, much less has any other computing system been subjected to so rigorous 
a test. 
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BIBLIOGRAPHY Z—XX]I 


9. Anon., “The Radar Range Calculator,” Bell Lab. Record, v. 29, 1951, 
p. 312-313. 


A circular slide rule calculator for computing the function of seven 
parameters which describes the effective range of a radar system. This 
function is a product of fractional powers of the parameters. 

F. J. M. 


10. A. A. CurrigE, “The general purpose analog computer,” Bell Lab. 
Record, v. 29, March 1951, p. 101-108. 


This article describes an analog computer using a basic three stage 
feedback amplifier for integration, differentiating and forming linear com- 
binations and with servo driven potentiometers for multiplication and 
finding functions like sin x. The computer contains 30 adders, 1 differen- 
tiator, 12 integrators, 9 servos which can drive 45 linear potentiometers 
and 15 function potentiometers, with 41 additional potentiometers which 
can be set by hand. There are six output recorders, two of which may be 
modified to be used as function input tables from graphs. The machine is 
set up by means of a three panel patch bay. 

F. J. M. 


11. L. L. Granpi & D. LEBELL, “‘Anaiogue computers solve complex prob- 
lems,”’ Radio and Television News, Radio Electronic Eng. ed., v. 46, 
No. 5, Nov. 1951, p. 70-71, 138-139. 


The Engineering Department of the Univ. of Calif. at Los Angeles has 
four continuous computers, an electrical differential analyzer, a mechanical 
differential analyzer, a network analyzer and a thermal analyzer. This 
article describes the problems for which each may be used and contains 
eight photographs including a color photograph on the cover. 

F. J. M. 


12. W. G. JAMEs, Logarithms in Instrumentation. United States Atomic 
Energy Commission ORNL-413, Oak Ridge National Laboratory, Oak 
Ridge, Tenn., 54 p. 

The author surveys the methods used for taking logarithms in instru- 
ments and indicates the accuracy and ranges available from the different 
methods. The logarithmic devices described include diodes with negative 
anodes, contact rectifiers, variable » tubes and potentiometers either wound 
with a special taper or with shunts between taps. 

F. J. M. 


13. L. H. Jonnson, Nomography and Empirical Equations. New York, 
Wiley; London, Chapman & Hall, 1952, ix + 150 p., 22.9 K 14.3 cm. 
Price $3.75. 


This book is an elementary text with the minimum of theory, and judged 
as such it is excellent. In particular the author considers questions of 
arranging nomograms for maximum accuracy. He also has a number of 
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good points in curve fitting, though his treatment of least square fitting 
could be improved a bit. 


Bell Laboratories 
Murray Hill, New Jersey 


R. HAMMING 


14. E. Lakatos, ‘Problem solving with the analog computer,” Bell Lab. 
Record, v. 29, March 1951, p. 109-114. 


A process is given for translating a differential equation into a form to 
which the differential analyzer is applicable. Two examples are briefly 
di d. 

eT F. J. M. 


15. Y. W. LEE & J. B. WiEsNER, “Correlation functions and communica- 
tion applications,’’ Electronics, v. 23, June 1950, p. 86-92. 


This article describes, largely in quantitative terms, some of the tech- 
niques and applications arising out of recent developments in the theory of 
communication, based on the statistical concept of information. The mean- 
ings of autocorrelation functions and cross-correlation functions are dis- 
cussed for stationary random processes, and these functions are applied to 
periodic processes. An electronic correlator, built at MIT, for obtaining and 
recording the autocorrelation curve of a random function is briefly described. 
The correlator obtains two sets of sampled values of the random function 
at regular intervals and forms pulses whose height and width are propor- 
tional to the two sets of sampled values, respectively. The autocorrelation 
function is found (approximately) by integrating the pulse areas. 

The problem of detection of periodic signals is discussed, and experi- 
mental results obtained with the correlator, showing a 30 db. gain for a 
signal which is 15 db. below the noise level, are presented. Additional appli- 
cations of the correlator to the determination of the impulse response of 
networks and the measurement of crosstalk in multichannel communication 
systems are indicated briefly. Finally, the importance of correlation func- 
tions in Wiener’s statistical prediction and filter theory is discussed. More 
quantitative information on many of the topics included in this article may 
be found in another paper by the same authors.' Details concerning the 
correlator are given in a technical report.? 


R. J. SCHWARZ 

Columbia University 
New York, New York 

1'Y. W. Lez, T. P. Cueatuam, JR., & J. B. WiEsNER, “Application of correlation analy- 
sis to the detection of periodic signals i in noise,” I. R. E., Proc., v. 38, 1950, p. 1165-1171. 

P. CHEATHAM, “Experimental determination of correlation functions and their 

ia to the statistical theory of communications,” Technical Report No. 122, Re- 
search Laboratory of Electronics, MIT. 


16. A. B. MAcNEE, “Some limitations on the accuracy of electronic differ- 
ential analyzers,” I. R. E., Proc., v. 40, 1952, p. 303-308. 


The author considers the errors that arise in solving a linear differential 
equation 


(1) > Ay” = F(t) 
n=0 
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with constant coefficients by means of an electronic analogue computer. 
Only adders and integrators are necessary. The steady state response of 
the simplest physically realizable adder is of the form 


Eout = K(z, + E; +--+ + E,)/(1 + juts). 


Similarly, for an integrator 


Eout/Ein = (1/jok JL1/(1 + 1/jouk) J[1/(1 + joT)] 


where u is the gain of the base amplifier and the last factor is due to the 
imperfect high frequency response of a physical integrator. Using these 
approximations, the author derives simple approximate expressions for the 
difference ¢,, ” = 1, ---, m of the characteristic roots s, (assumed distinct) 
of (1) and the characteristic roots s,’ of the differential analyzer. The 
m+ 1 additional roots of the machine solution [introduced by the 
“1/(1 + jwT)” factors which raise the order of (1) ] are, under physically 
reasonable assumptions, widely separated from the s, roots and have nega- 
tive real parts. Numerical examples are included. The case in which (1) has 
multiple roots is briefly touched upon; the phenomena which arise when 
the Jordan normal form is not diagonal are not mentioned. 


K. S. MILLER 
New York University 


New York, New York 


17. R. E. Scott, An analog device for solving the approximation problem of 
network synthesis. Technical Report No. 137, June 8, 1950, Research 
Laboratory of Electronics, Massachusetts Institute of Technology, Cam- 
bridge, Mass. 


Mathematically the problem treated here is that in which one endeavors 
to find a rational function which is a good approximation or the imaginary 
axis to a prescribed complex valued function. The desired function is the 
‘network response” function from which a network can be synthesized with 
a given frequency response corresponding to the prescribed values along 
the imaginary axis. The rational function must satisfy a number of condi- 
tions, for instance, the numerator and denominator must be expressible 
as real polynomials and the roots of the denominator must all be in the 
left half of the complex plane. It is possible to represent the logarithm of a 
rational function with a prescribed set of zeros and poles as the potential 
function of a planar flow of electric current with corresponding sources 
and sinks. 

The device described in this report has as its conducting plane a sheet 
of ‘‘Teledeltos” paper. Current sources are provided and also probes which 
are stuck into the paper to represent the sources and sinks. The values of 
the potential function along the imaginary axis are sampled from fixed 
contacts by a commutator and displayed on an oscilloscope. The movable 
probes are adjusted until the desired behavior is obtained. One major 
problem in this process is due to the finite extent of the conductor. A num- 
ber of solutions for this have been proposed.' In this device, a logarithmic 
transformation is made which with certain symmetry considerations permits 
one to substitute a strip instead of the infinite plane. The strip is termi- 
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nated at each end at some distance from the computing portion to represent 
zero and infinity. To obtain the equivalent of phase the fixed contacts along 
the imaginary axis are doubled so that potential slope can be indicated. 
The report contains considerable information on both the construction and 
use of this device and an analysis of the effect of individual errors. A photo- 
graph shows the entire machine mounted on a laboratory cart. 

F. J. M. 


1Cf. A. R. BoorHroyp, E. C. CaHerry & R. Makar, Inst. Elect. Eng., Proc., v. 96, 
1949, p. 163-177; MTAC, v. 5, p. 49-50. 


18. J. R. SHan & L. Jacoss, “Investigation of field distributions in sym- 
metrical electron lens,” Jn. Appl. Physics, v. 22, 1951, p. 1236-1241. 


The potential field for the electron lens is obtained by the use of an 
electrolytic tank. It was determined experimentally that an impedance of 
30 megohms for the probe was necessary. With these methods, the crossing 
of the nodes at the central symmetry point was accurate to within about a 
quarter of a degree. The field obtained from the electrolytic tank was com- 
pared with that obtained by a relaxation method. The maximum difference 
is about 1.5%. 

F. J. M. 


19. B. A. Soxo.orr, “Principe et réalisation d’une machine mathématique 
dite ‘Operateur Mathématique Electronique’ OME,” Annales des Télé- 
communications, v. 5, no. 4, April 1950, p. 143-159. 


The mathematical machine described is an electronic differential ana- 
lyzer with the customary feedback amplifier integration and addition. The 
output indicator is an oscilloscope and camera combination. Multiplication 
is by servo driven potentiometers. A discussion of stability is given for the 
constant coefficient case of linear differential equations, including the use 
of the determinantal conditions for positive definiteness. 


F. J. M. 


20. G. J. TaAuxE & R. L. Stoker, “Analytical studies in the suppression of 
wood fires,’”” A. S. M. E., Trans., v. 73, 1951, p. 1005-1020. 


Different methods of suppression of wood fires are analyzed thermally 
using electrical analogy methods. The latter effect the solution of the heat 
equation by means of an R-C network. The method has been previously 
described. ! 

JEROME H. WEINER 
Columbia University 
New York 27, New York 


1V. Pascuxis & H. D. Baker, “A method for determining unsteady-state heat 
transfer by means of an electrical analogy,’ A. S. M. E., Trans., v. 64, 1942, p. 105-110. 


NOTES 


139. THE RAND COLLECTION OF ILLUSTRATIVE APPROXIMATIONS.— 
In recent years, the Numerical Analysis Department at RAND has been 
preparing loose leaf sheets that contain interesting and useful approxima- 











252 NOTES 


tions to a number of the higher transcendental functions. The data from a 
sample sheet are displayed at the end of this note together with a photo- 
graphic reproduction of the error curve appearing on this sheet. Some sixty 
odd sheets in this series have now been prepared, and up to date collec- 
tions of these sheets have been distributed to some three hundred people 
working in the field of numerical analysis throughout the United States 
and in a few foreign countries. 

The approximations given in this series of loose leaf sheets are of both a 
practical and of an illustrative nature. In a high speed digital machine, 
these approximations may take the place of bulky tables or the place of 
awkward series developments. For the hand computer, an easily evaluated 
expression will often be of use when a required table is unavailable for 
consultation. Beyond this, however, the sheets will be useful for the insight 
they give the practical computer in the approximation of functions. Starting 
with approximations concerning the common logarithm, sheets for an ever 
increasing number of the useful transcendental functions have been pre- 
pared and distributed. Sheets of approximations in this series are already 
available concerning the logarithmic function, the exponential function, the 
inverse tangent, the sine function, the inverse sine, the Gamma function, 
the Gaussian error integral, the inverse Gaussian error integral, the com- 
plete elliptic integrals of first and second kind, the exponential integral 
and a number of other special functions. 

As the approximations in this series are intended to be illustrative as 
well as practical, great care has been taken in the accurate leveling of the 
error curves. The sheets of approximations thus give interesting information 
concerning the location of roots—points at which function and approxima- 
tion agree—and the location of extremals—points at which the deviation 
between function and approximation is locally a maximum in an absolute 
or relative sense. The accurate and carefully drawn error curves that appear 
on each sheet will do much to give the reader a feel for the nature of ana- 
lytical approximation in a wide variety of typical cases of practical impor- 
tance. In order to provide this important information, the coefficients in 
the approximations given have been recorded to perhaps two, three or four 
more decimals than would be required for practical considerations of utility. 
Should the coefficients given be of awkward size for use in a given computing 
machine, it is permissible to round the numbers quite severely. Of course 
the beauty and character of the error curve will be lost as a result of such 
mistreatment, but this is not a matter of practical concern in the use of the 
approximation. 


CEcIL HASTINGS, JR. 
The RAND Corporation 


Santa Monica, California 


The RAND Corporation 
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140. AN ALTERNATIVE ‘‘PREDICTOR-CORRECTOR” PROcEss.—The pur- 
pose of this note is to report the successful use, on a system of 3 non-linear 
ordinary differential equations of order 2 (with initial conditions), of a pair 
of formulas analogous to Milne’s.! 


The system of differential equations is of the form 
(1) £; = fi(di, Le, Ls, X1, X2, Xs, 4), (4 = 1, 2, 3). 
The formulas used for prediction are 
(2) 2) = 54,0-) — 47; 4 2h[4;°-” + 22; ], (¢ = 1, 2, 3). 


The correction formulas are 

(3) &D = oe 4 t [Sé:+ 4+ 8#.~ —4,-D], (i = 1, 2, 3) 
and their companions 

(4) x) = x 4 4 [Se 4+ 82, — 4-D], (¢ = 1,2, 3). 


In the above, h denotes the integration step. 
The process is as follows: Given two sets of ‘‘starting’’ values, x;“~", 
:-, €;°-) and x;™, 2; #;™ [the values at ¢ = tp. andt, = tnt h, 
resp. ], the first approximations ,¢;+" are predicted by (2), then (4) are 
used to obtain ,x;“"+!), These first approximations are substituted in (1) to 
obtain ,@;+”, which are then used in (3) to obtain »¢;*+”. The above 
steps may be repeated, to give 3¢;"+, etc. If the @;*" prove acceptable 
(see below), then ;x;“*+” are computed and accepted. 
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An estimate of the error in the second approximation to #;**” is avail- 
able from the error terms of (2) and (3), viz. h4x,;°(£)/6 and —h*x;*(n)/24, 
respectively, where both ~ and 7 are on the interval tn; < ¢ < tay: and 
where x;* denotes the fifth derivative of x;. If we let Diy = justi) — ja», 
then, as in Milne,' we find that the error in »¢;*" is approximately Dj,/5. 
Thus, if all the |Dji/5| are insignificant, one iteration is probably sufficient, 
whereas if any |Dj/5| is intolerably large, a smaller h and/or more itera- 
tions are called for. 

In the problem to which the process was applied (it was carried out on 
our IBM CPC), h = 0.1 was chosen initially. This proved to be a good 
choice, with 2 iterations, from ¢ = 0 to t = 4.0, as shown by Dy and Dj, 
which were all listed for monitoring. In most instances, some |Dj/5| was 
too large (>4-10-*), whereas no |Dj2| exceeded 10-*, thus indicating that 
further iterations would probably not affect the 3rd decimal place. The 
combination h = 0.1, and two iterations, was actually used to ¢t = 4.4, at 
which place some | Dj2| increased significantly ; the # was then decreased to 
0.05 and the process continued from ¢ = 4.0 with two iterations (although 
the previously computed values at ¢ = 4.3 were probably satisfactory, it 
was felt best, for the sake of smoothness and safety, to back up to ¢ = 4.0 
before changing to h = 0.05). The new combination gave good results as 
far as t = 5.15, when some |Dj2| again became significant. Here it was 
decided to use h = 0.05 again, but increase to 4 iterations; this worked 
successfully as far ast = 4.90. 

It may be of interest to note that the need for smaller 4 and/or more 
iterations arose because the f; in (1) each would (apparently) have become 
infinite in the neighborhood of t = 4.7. The results were checked by differ- 
encing the (accepted) values of the x;. In addition, there was an internal 
check available in the problem, viz. that z,? + a? + #32 = K. Interestingly 
enough, this check failed approximately at the same time that the last 
computed |D,;| became intolerably large. 

One of the advantages of this process is that it requires only two sets of 
“starting” values rather than four. This means that in case it becomes 
necessary to decrease h at some point of the process, only one additional 
“point” is needed. This added point can be calculated, e.g., by means of 
the modified Euler method,? as can the additional point needed initially. 
Another “‘advantage”’ (rather a nebulous one, actually, since it is based on 
esthetics rather than mathematics) is that one has a feeling of greater con- 
fidence in values based on others in their more immediate neighborhood 
(two starting values as opposed to four). Obviously, examples exist which 
would vitiate this belief. 

T. H. SouTHARD 


E. C. YOWELL 
National Bureau of Standards 
Los Angeles, 24 


The preparation of this paper was sponsored (in part) by the Flight Research Labora- 
tory, Wright Air Development Center, USAF. 


1W. E. MILNE, Numerical Calculus. Princeton, 1949, p. 135. 
2 J. B. ScarsorouGH, Numerical Mathematical Analysis. Baltimore, 1950, p. 235. 
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141. A NEw CALCULATION OF EULER’s CoNSTANT.—In 1878 J. C. ADAMS 
published his notable paper! on the evaluation of Euler’s constant. By two 
separate calculations he determined approximations which agreed to 263D. 

Since that time no attempt seems to have been made to extend the 
approximation beyond that point, until a recent study of Adams’ work 
induced me to undertake the task. 

Following the procedure adopted by Adams, I used the Euler-Maclaurin 
summation formula applied to the harmonic series, which yields the equation 


n 1 m (- 1)"1B, 
_— me. oe an Pe ef 
y=Lk ins an + 2rn* 


k=1 r=1 


+ Rau. 


For fixed n, the minimum absolute value of R,, is nearly equal to the numeri- 
cal value of the term of the asymptotic series given by taking m = [nx — 3]. 
For this calculation of y it was decided to take m = 1000, principally 
because of the resultant ease of evaluating the sum of the asymptotic series 
to about 350D. 
The evaluation of the sum of the first thousand terms of the harmonic 
series to 350D was accomplished by means of an artifice used by Adams. 


168 
He replaced that sum by an equivalent, namely —43 + > a:p;-*, where 
i=1 


p® < 1000 < p*, k; > 1, and a; is a positive integer <p/* determined 
uniquely by decomposing the terms of the partial sum of the harmonic 
series into partial fractions with prime-power denominators, adding all the 
component fractions involving the same prime ;, and finally reducing the 
numerator modulo p;*. This part of the calculation also was facilitated by 
the use of an accurate manuscript table of the complete decimal periods of 
the reciprocals of all primes less than 1000, which I had previously calcu- 
lated in order to check Gauss’ table of circulating decimals.” 

The only remaining datum required was In 1000. This was deduced 
from the 330D approximation to In 10 calculated by H. S. UHLER.* 

Inasmuch as Uhler guaranteed his value of In 10 to within 2 units in 
the 329th decimal place, and all other parts of this calculation were care- 
fully checked, it is believed that the following approximation to Euler’s 
constant is correct to at least 328D. It confirms the accuracy of Adams’ 
approximation to 262D. 


y =0.57721 56649 01532 86060 65120 90082 40243 10421 59335 93992 
35988 05767 23488 48677 26777 66467 09369 47063 29174 67495 
14631 44724 98070 82480 96050 40144 86542 83622 41739 97644 
92353 62535 00333 74293 73377 37673 94279 25952 58247 09491 
60087 35203 94816 56708 53233 15177 66115 28621 19950 15079 
84793 74508 57057 40029 92135 47861 46694 02960 43254 21519 
05877 55352 67331 39925 40129 674(28) 

J. W. WRENCH, JR. 

4711 Davenport St., N.W. 

Washington 16, D. C. 


1 Roy. Soc. London, Proc., v. 27, 1878, p. 88-94. 

References to earlier calculations are given in FMR, Index, p. 95. 
2 MTAC, v. 4, 1950, p. 222, 223. 
3 Nat. Acad. Sci., Proc., v. 26, 1940, p. 205-212. 











256 NOTES 


142. THE DETERMINATION OF A LARGE PRIME. [EDITORIAL NOTE: The 
primality of (2! + 1)/17 as established by A. FERRIER was very briefly 
announced in MTAC, v. 6, p. 61. Since this is probably the last ‘‘largest”’ 
prime to be identified by hand computing methods (primes like 2!7° — 1 
would require more than a century of desk calculator work to establish by 
any known method), it may be of interest to give the following details 
quoted from Ferrier’s letter of July 14, 1951. ] 

I have established that 


N = (28 + 1)/17 = 2098 89366 57440 58648 61512 
64256 61022 25938 63921 


is a prime and is in fact the largest prime known. The method used is the 
following. 
I first found that NV divides 3" — 3". In fact for m = 2'* we have 


3” = —566 68397 79443 57425 64352 
U (mod N) 


Il 


and 
U? 


] 


32n 


ll 


—9 (mod N). 
Hence 
3-(3")6 = 3N = 3-(—9)8 = 3” (mod JN). 


That is, N divides 3" — 3". Next we see that 
N —1 = (28 — 24)/17 = 24(27 — 1)(27 + 1)/17. 
The largest prime factor of N — 1 is that dividing 2” + 1 and is 
pb = 48 78248 87233. 


Writing VN — 1 = pm we find that 3’ — 1 is prime to N. Applying the 
theorem quoted and corrected in MTAC, v. 3, p. 497, and v. 5, p. 259, 
respectively, it follows that every prime factor of N is of the form px + 1 
as well as 296x + 1, that is of the combined form 


14439 61666 20968x + 1 = gx + 1. 


For g = 1(1)11, gx + 1 is divisible by small primes. Hence every prime 
factor of N exceeds 12g. 
Writing N = A? — B? we have 


2A < 12g + N/(12q) < 1.3-10*. 
However 
2A N + 1 (mod q’) 


1885 97808 71263 54966 31614 94450 (mod g?) 


so that 2A > 1.8-10%. 
Thus A does not exist and JN is a prime. 


Ii 


A. FERRIER 
Collége de Cusset 
Allier, France 
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QUERIES—REPLIES 


50. A DEFINITE INTEGRAL (Q 41, v. 6, p. 125). 
In the definite integral 


[= f u— exp (—2u — u~*)du 
0 
put z = 2x*, v = xu so that 
a 
[= f v exp [—x?(2v + v~*) ]dv. 
0 
In this form the integral can be evaluated asymptotically by LAPLACE’s 


method.' f(v) = 20 + v~? has its maximum at v = 1, and for large x the 
integral is approximated by 


2r ; 
lara | exp [—x?f(1) ] 


(3/3) ix-te-**, 


or 


This approximation may be used to start an asymptotic expansion: suc- 
cessive terms may be computed from the differential equation stated in the 
query. 


A. E. 
1See D. V. WippER, The Laplace Transform. Princeton, 1946, p. 277-280. 


EpitoriAL Note: Dr. J. ERNEsT WILKINS also obtained the dominant 
term of the asymptotic expansion of this integral for large values of 2. 


. 7, for 1950 read 1949. 

. 9, for 21 read 206. 

. —4, for yi read 1. 

. 18, for 33 read 5. 

58, |. 20, for 35 read 5. 

61, 1. —8, for 116 read 118. 

82, 1. 6, for TaBListy read TaBLifSy. 
. 86, 1. —17, for p = 0.6 read p = 0.0. 
161, 1. —11, for 7956 read 7556. 


CORRIGENDA 
3, p. 355, 1. —12, for 10% read 10%. 
6, p. 20, 1. —5, for n < 50 read n < 50. 
6, p. 25, 1. 14, for 1.2(.5) read 1.2(.05). 
6, p. 25, 1. 14, for 2-3D read 3-4D. 
6, p. 25, 1. —20, for R52 read R53. 
6, p. 32, 1. —7, for JOHNSON read JOHNSTON. 
l 
I 


ADAAAAAAG 
CeCrtrerrres % 


aun ww 
ant dw 


€$#<¢<¢<¢e4eee6464868 48446484 46486 








SUBJECT INDEX 


ABC 185 

Aberdeen Proving Ground Computers 50, 
113 

ACE 186 

Aerodynamics see Fluid Mechanics 

Airy Integral 98 

Algebraic Equations 21, 82-83 

Approximating Functions 251-253 

Arc Length 81 

Arc Tangent 19 

Arc Tangent Relations 19 

Astronomy 162 

Astrophysics 230-231 

Automatic Computing Engine 186 


Bell Telephone Computers 183 

Benedict Equation 114 

Bernoulli Numbers 151 

Bessel Functions 24-30, 39, 95, 99-100, 159, 
237 

Bessel Functions, Incomplete 41 

Bessel Functions, Integrals of 159 

Bessel Functions, Roots of 24-25, 28, 32, 
58-59, 97-98 

Beta Function, Incomplete 36 

Binary Automatic Computer (BINAC) 183 

Binomial Coefficients 18 

Binomial Distribution 24 

Biography 125 

Blasius Equation 233 

Boolean Algebra 243-244 

Burroughs Laboratory Computer 113, 117 


CADAC 113, 186 

Calculating Machines see Computing Ma- 
chines 

Chi Square Distribution 21-22, 36, 56 

Circle Computer 50, 186 

Coding 49-52, 112-113, 167-182, 186 

Commonwealth Scientific and Industrial 
Research Organization Computers 53, 
167-172 

Complex Numbers 167 

Computing Machines, Analogue 1-8, 55-58, 
119-123, 189-204, 241, 246, 251 

Computing Machines, Digital 41-54, 102- 
119, 167-190, 238-247 

Congruences 167 

Consolidated Electronic Digital Computer 
186 

Constants 19, 31-32, 151 

Continued Fractions 127-133 
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Conversion Equipment 49, 54, 114 

Conversion Tables 18-19, 167 

Correlation 85-86, 229, 249 

Cosine, Hyperbolic 151 

Cosine Integral 151 

Cosines, Cyclotomic 102 

Cotangent 163-166 

Coulomb Functions 92-95 

CPC 114, 167 

CRC 105 Decimal Digital Differential Ana- 
lyzer (DDDA) 111-112, 247 

CSIRO, Mark I, 53, 167-172 

Cubic Equations 82-83 

Cubic Fields 21 

Cylinder Function 97 


Data Processing 114, 186 

Decimal Equivalent of Fractions 101 

Decimals, Circulating 255 

Determinants 68-72 

Difference Analyzer 1-8 

Differences 212-214 

Differential Analyzers 41-49, 53, 56-57, 
102-112, 246-247, 249-250 

Differential Equation Solvers 1-8, 55-56 

Differential Equations 1, 4-5, 29-30, 162 

Differential Equations, Numerical Solution 
of 63-78, 138-150, 152-153, 219-223, 234- 
236, 253-254 

Digital Differential Analyzers 41-49, 102- 
112, 247 

Diophantine Equations 20-21, 245 

Dirichlet’s Singular Integral 28 

Distribution Functions 23-24, 40, 84-91, 
154-157, 227-229 

Distribution of Digits 19, 124 


e 32, 123-124 

EDSAC 49-52, 61, 118 

EDVAC 113, 115 

Eigenvalues 11-16, 26, 29, 54, 196-199 

Electromagnetic Waves 94 

Electronic Computer Corporation Comput- 
ers 50 

Electronic Delay Storage Automatic Cal- 
culator (EDSAC) 49-52, 61, 118 

Electronic Discrete Variable Computer 
(EDVAC) 113, 115 

Electronic Numerical Integrator and Com- 
puter (ENIAC) 32, 95, 113, 183, 239 

Elecom 100 Computer 187 

Elliptic Functions 159, 236 
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Elliptic Integrals 28, 160 

Engineering Research Associates Ccm- 
puter 118, 172-183 

ENIAC 32, 95, 113, 183, 239 

Equation Solvers 54, 119, 121-122 

Equations, Algebraic see Algebraic Equa- 
tions 

Equations, Linear see Linear Equations 

Equations, Transcendental see Transcen- 
dental Equations 

ERA 1101 Computer 118, 172-183 

Errata 21, 24, 35-39, 99-102, 151, 159, 163- 
166, 227, 231, 236-238 

Error Function see Probability Integral 

Euler’s Constant 255 

Euler-Maclaurin Summation Formula 255 

Evolutes 224 

Exponential Function 84, 159 

Exponential Integral 28, 94, 151, 236 


Factor Tables 39 

Factorials 32-33, 151 

Factorizations 39, 59-60, 102 

Ferranti Computer 50 

Fibonacci Numbers 20 

Fields, Cubic see Cubic Fields 

Fluid Mechanics 30-31, 33-34, 41, 93-94, 
160-161, 233 

Fourier Coefficients 26 

Fourier Series 29-30, 205 

Fourier Synthesizer 123 

Fresnel Integrals 163 

Functions (special) see Bessel, Beta, Cou- 
lomb, Cylinder, Distribution, Elliptic, 
Error, Exponential, Gamma, Hankel, 
Hyperbolic, Hypergeometric, Hyperel- 
liptic, Legendre, Logarithmic, Mathieu, 
Partition, Probability, Riemann Zeta, 
Struve, Symmetric, Trigonometric, 
Wave, Zeta 


Gamma Function 28, 32-33, 151, 159, 252 
Gamma Function, Incomplete 21, 36 
Gauss Quadrature 225 

Gauss Transformation 9-10, 14-15 

Gear Ratios 101 

Goldbach Problem 238 

Groups 225 


Hankel Functions 28, 97, 237 
Harmonic Analyzer 56, 123 
Harmonic Synthesis 123 
Harvard Mark III 118, 183 


Helmholtz Equation 29 

Hermite Interpolation Formula 225-227 

Hermite Polynomials 225-227, 232 

Hyperbolic Functions 28-29, 151, 159 

Hyperbolic Functions, Inverse 159 

Hyperelliptic Functions 93 

Hypergeometric Functions 30-31, 159-160 

Hypergeometric Functions, Confluent 94- 
95, 155, 159-161 

Hypergeometric Integrals 28 


Inequalities 31-32 

Institute for Advanced Study Computer 
245-246 

Integral Equations 133-138, 230-231 

Integrals, Definite 24, 26-28, 40-41, 92, 125, 
231, 233, 235-237, 257 

Integrals (special) see Cosine, Elliptic, 
Exponential, Fresnel, Probability, Sine 

Integrators 123 

International Business Machines 40, 114— 
118, 120, 124-125, 182, 184, 190, 254 

Interpolation 212-214 

Interpolation Coefficients 225-227 


Jacobs Instrument Company Computers 
186, 246 


Kloosterman Sums 102 


Laguerre Polynomiais 162 

Laplace Transforms 28, 100-101, 158-159 
Least Squares 39-40, 228 

Legendre Functions 30, 159 

Legendre Functions, Associated 30, 98-99 
Legendre Polynomials 95-97 

Linear Equations 9-17, 34, 54-55, 190-203 
Linear Equation Solvers 55 

Logarithmic Function 159 

Logarithms 150-151 


Magnetic Drum Digital Differential Ana- 
lyzer (MADDIDA) 49, 114 

Manchester University Digital Computer 
241-242 

Mathieu Functions 25-27, 29-30, 159 

Matrices 10-17, 34, 54, 78-81, 116 

Mechanical Quadrature see Numerical In- 
tegration 

Memory Devices 243 

Ménage Numbers 152 

Mersenne Primes 61, 125, 205 

Mock Theta 224 
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Monte Carlo Method 78-81, 115, 138 

Moore School Automatic 
(MSAC) 113, 186 

Multiplication Tables 61, 205 


Computer 


National Bureau of Standards Eastern Au- 
tomatic Computer (SEAC) 50, 113, 
118, 183-185 

National Bureau of Standards Western 
Automatic Computer (SWAC) 50, 61, 
113, 186, 189, 205 

Naval Proving Ground Calculators 50, 113 

Network Analyzer 56-57, 250-251 

Nomograms 123, 248-249 

Non-Linear Equations 114-116 

Number Theory 19-21, 39-40, 59-60, 81-82, 
151-152, 224, 238, 245, 256 

Numerical Integrations 63-78, 
215-219, 220-221, 225-227 

Numerical Solution of Differential Equa- 
tions 63-78, 138-150, 152-153, 219-223, 
234-236, 253-254 


138-150, 


OARAC 113 

Office of Naval Research Relay Computer 
185, 207-212 

Operational Calculus 100-101, 158-159, 162 

ORDVAC 50, 113, 118 

Orthogonal Polynomials 28, 40 

Output Devices 244-245 


Parabolic Coordinates 81 

Partition Function 20, 151-152 

Perfect Numbers 61 

Physical Tables 56 

Pi 101 

Planimeters 120 

Poisson Distribution 84 

Polar Coordinates 224 

Polynomials see Hermite, Laguerre, Or- 
thogonal 

Primes, Large 61, 256 

Primes, Lists of 38, 40, 81-82, 163, 167 

Probability Functions see Distribution 
Functions 

Probability Integral 28, 94, 102, 151, 159, 
232, 252 

Probability Integral, Derivatives of 232 

Programming 49-52, 112-113, 167-182, 186 

Punch Card Methods 56, 114-117, 133-138, 
182 

Punch Card Tables 167, 204-205 


Quadratic Arc Computer (QUAC) 187 
Quintic Equations 21 


Radiation 29-30 

Random Walks 78-81 

Raytheon Computer 50 

REAC 190-203, 246-247 

Reciprocals of Complex Numbers 167 

Reciprocals of Factorials 32 

Reeves Electronic Analogue 
(REAC) 190-203, 246-247 

Regressions 87-88, 156 

Relaxation Methods 9-11 

Remington Rand Computers 187 

Riccati Bessel Functions 95-96 

Riemann Zeta Function 161-162 

Rounding Off Errors 212-214 


Computer 


Scattering Function 96 

Schroedinger Equation 29 

SEAC 50, 113, 118, 183-185 

Selective Sequence Electronic Computer 
(SSEC) 183-184 

Simpson’s Rule 222-223 

Sine, Hyperbolic 151 

Sine Integral 151, 237 

Slide Rule 56, 121 

Spherical Harmonics 29 

Square Roots 19, 57 

SSEC 183-184 

Statistics 21-24, 35-38, 83-91, 153-158, 163, 
227-230 

Struve Functions 237 

Student’s ¢ 89 

Surveying 151 

SWAC 50, 61, 113, 186, 189, 205 

Symmetric Functions 224-225 

Symmetric Group 225 


Tangent 163-166 

Tangent, Inverse see Arc Tangent 

Taylor’s Series 218 

Telegraph Equation 231-232 

Thermodynamics 114-116 

Theta Functions 159 

Transcendental Equations 32-33, 204 

Trigonometric Functions 26, 29, 32, 100, 
151, 159, 163-166, 224 

Trigonometric Functions, Inverse 159 

Universal Automatic Computer (UNI- 
VAC) 117, 119, 183, 186, 247 

University of Manchester Maciiine 118 
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University of Toronto Electronic Com- 
puter (UTEC) 50 


Wave Functions 29, 231, 233-234 


Whirlwind I, 113, 118, 185 
Wilson’s Quotient 238 


Zeta Function 161-162 
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